Entanglement renormalization and gauge symmetry 



O 

Oh 
< 



L. Tagliacozzo and G. Vidal 
School of Mathematics and Physics, the University of Queensland, QLD 4072, Australia 

(Dated: April 12, 2011) 

A lattice gauge theory is described by a redundantly large vector space that is subject to local 
constraints, and can be regarded as the low energy limit of an extended lattice model with a local 
symmetry. We propose a numerical coarse-graining scheme to produce low energy, effective descrip- 
tions of lattice models with a local symmetry, such that the local symmetry is exactly preserved 
during coarse-graining. Our approach results in a variational ansatz for the ground state(s) and 
low energy excitations of such models and, by extension, of lattice gauge theories. This ansatz 
incorporates the local symmetry in its structure, and exploits it to obtain a significant reduction of 
computational costs. We test the approach in the context of a Z^ lattice gauge theory formulated 
as the low energy theory of a specific regime of the toric code with a magnetic field, for lattices with 
up to 16 x 16 sites (16 2 x 2 = 512 spins) on a torus. We reproduce the well-known ground state 
phase diagram of the model, consisting of a deconfined and spin polarized phases separated by a 
continuous quantum phase transition, and obtain accurate estimates of energy gaps, ground state 
fidelities, Wilson loops, and several other quantities. 

PACS numbers: 03.67.-a, 05.50.+q, 11.25.Hf 
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I. INTRODUCTION 

Gauge invariance is one of the most important concepts 
in modern physics. It is at the core of the equivalence 
principle of general relativitjff^ as well as an essential 
ingredient in the quantum field formulation of the stan- 
dard model of particle physical] More broadly, gauge 
theories are useful in several research areas, including 
condensed matter, nuclear and high energy physics. 

Wilson's lattice gauge theory-^ was proposed in 
order to study quantum chromodynamics QCD non- 
perturbatively. By replacing continuous space-time by 
a discrete lattice, it offers both a specific regulariza- 
tion scheme and a convenient starting point for numer- 
ical studies. Through large scale Monte Carlo compu- 
tations, lattice gauge theory has provided, amo ng o ther 
results, numerical evidence of QCD confmementP^ and 
of the presence of a chiral condensate^^, the determi- 
nation from first principles of the value of the quark 
massesD^ and numerical evidences of the for mation of 
quark-gluon plasma at high temperatures 119 . In spite 
of their incontestable importance in lattice gauge theory 
calculations, Monte Carlo sampling techniques can not 
be applied to systems with large chemical potential due 
to the fermionic sign problem. It is therefore desirable to 
explore alternative approaches. 

In the last few years, tensor network algorithms 
have received increasing attention as variational, non- 
perturbative methods to study spin lattice systems. 
The simplest tensor networ k variational ansatz, the 
matrix product state^^H (MPS), is the basis of 
White's extr emely s uccessful density matrix renormal- 
ization g r0U rji23ftl25| (DMRG) algorithm to compute the 
ground state of spin chains. Several generalizations of 
the MPS to two and larger spatial dimensions exist 
both for spin and fermionic s yste ms, such as the pro- 
jected entangled-pair state9^5HI32|(PEPS), a i so referred 



to as tensor product states, a nd the multi-scale entan- 
glement renormalization ansatzSHHUni (MERA), both of 
which lead to scalable simulations. Since tensor net- 
work algorithms do not suffer from the sign problem, 
they are suitable to study problems beyond the reach 
of Monte Carlo techniques. Recent ca lculations in- 
clude frustr ated an tiferromagnet d 131 * 147 ^^ and interact- 
ing fermion d 136 ! 151 ! in two spatial dimensions. 

The goal of this paper is to explore the use of tensor 
network techniques within the Hamiltonian formulation 
of lattice gauge theory. We will consider the simplest 
non-trivial case, n amely Zi lattice gauge theory in two 
spatial dimension d 1 1 5 I 1 52 H1MI ; w Jth a Hamiltonian that 
contains the usual kinetic and potential terms for the 
gauge field and no fermionic matter. The ground state 
phase diagram of the theory, which we aim to reproduce, 
is already well understood, due e.g. to previous studies 
using Monte Carlo techniques. Therefore the aim of the 
present paper is not to uncover new physics, but rather 
to describe a new approach and to confirm its validity in 
a well-known, simple context. The merit of the present 
strategy resides in that it can be generalized to more 
complex settings beyond the reach of Monte Carlo tech- 
niques. 

In this work we will describe (i) a coarse-graining 
scheme for the lattice model specifically designed to pre- 
serve and exploit its local Zi symmetry; (ii) a variational 
ansatz for the ground state(s) and low energy states of 
theory; and (iii) numerical results to illustrate the po- 
tential of the approach. We will regard lattice gauge 
theories from the broader and more recent perspective of 
topological order^H The Z^ lattice gauge theory can be 
studied as the low energ y sec tor of a spin model, namely 
Kitaev's toric code modeP^, in a suitable limit (see also 
Refs. 157 160). The hopping term of the gauge field is 
implemented by deforming the toric code with a mag- 
netic field, as investigated in Refs. 11611 and 11621 In the 



deformed toric code model, the Z 2 group is a local sym- 
metry of the spin Hamiltonian, and states are not forced 
to be gauge invariant. In particular, excited states exist 
that are gauge covariant. Such states can be understood 
as describing the presence of static matter in specific lo- 
cations of the lattice. Our variational ansatz can also 
represent such states. 



A. Previous work 

Previous related work can be div ided into two cate- 
gories. On the one hand, Sugihara^J has investigated the 
use of a MPS to describe the ground state of Z 2 lattice 
gauge theory in quasi-one dimensional systems, both di- 
rectly or through a mapping to the quantum Ising model. 
We notice that no attempt is made to adapt the MPS, 
primarily a tensor network for one-dimensional systems 
(where Z 2 lattice gauge theory is trivial) to the presence 
of the symmetry 

On the other hand, Schuc h, C irac and Perez-GarciaPSl, 
as well as Swingle and WertlSE have recently considered 
a PEPS representation of a state of a two dimensional 
lattice with topological order and have investigated how 
to e xtra ct topological information from it, while Chen 
et alP^SI have stressed the importance of explicitly pre- 
serving certain symmetry in a PEPS in order to repre- 
sent topological order. In contrast with the present work, 
where we propose and benchmark an algorithm to com- 
pute the ground state(s) of a two dimensional system 
with topological order (in some limiting regime) or, re- 
latcdly, of Z 2 lattice gauge theory, none of these previous 
contributions explains how to actually obtain, given a 
Hamiltonian of a topologically ordered system, a tensor 
network representation of its ground state(s). 



B. Guide 




FIG. 1. Entanglement entropy of half a 4 x 4 (32 spin) lattice 
jC in the ground state |$+,+) of H^ c , as a function of the mag 
netic field h x . The critical magnetic field h c x Tit = 0.3285(1; 
is denoted by a vertical line. The blue (upper) curve is a lower 
bound to the entanglement entropy of the deformed toric code 
model without Wexact- For comparison, the green (middle) 
curve corresponds to the entanglement entropy of the ground 
state of the quantum Ising model on a 4 x 4 (16 spin) lattice 
with periodic boundary conditions and no vacancy. Finally, 
the red (lower) curve corresponds to the entanglement left in 
the ground state of H^ c after applying the analytical trans- 
formation H4xact- Notice that the amount of entanglement 
is very similar to that of the Ising model. Similar reduction 
of the entanglement entropy was obs erved b y several authors 
in the context of global symmetrie s ' 168 ! 16 ^ Thus, by imple- 
menting a local version of the duality transformation to the 
Ising model, the analytical transformation W eX act maps the 
robustly entangled ground state 1$+,+) of the deformed toric 
code model to a significantly less entangled ground state. 



A reader who simply wants to obtain an overall picture 
of the approach may prefer to read the summary below 
and then have a look at the benchmark results of Sect. 



ED 



The paper is organized in several sections. Sects. \U\ 
|IV| are entirely devoted to introducing background mate- 
rial, whereas Sect. |V| describes the proposed approach. 
Sect. VI presents benchmark results, and Sect. VII con- 



cludes with a discussion of the proposed technique and 
of a number of possible generalizations to be pursued in 
future work. 

In more detail, Sect, [n] reviews the toric code and 
its magnetic field deformation, Sect. |III| reviews two- 
dimensional Z 2 lattice gauge theory and its duality with 
the two-dimensional quantum Ising model, and Sect. |IV| 
contains a short introduction to coarse-graining transfor- 
mations based on entanglement renormalization and to 
the MERA. A reader with sufficient previous knowledge 
on the toric code, Z 2 lattice gauge theory and entangle- 
ment renormalization may decide to skip Sects. H]|IV in 
a first reading and jump directly to Sect. |Vj which de- 
scribes our proposal — possibly returning to the Sect. 
IllllIVI to seek clarification on nomenclature. 



U. summary 

As explained in Sect. [Vj the lattice model is coarse- 
grained by means of a transformation that breaks into 
two parts Wexact and W num . The first part Wexact can 
be viewed as a local version of the duality transforma- 
tion to the quantum Ising model. The spin model is 
divided into blocks of spins, and within each block a du- 
ality transformation to the Ising model is implemented. 
This produces two types of spins: at the boundary be- 
tween blocks, spins which have not yet been mapped to 
the Ising model (referred to as constrained spins) and in 
the interior of each block, spins that have been mapped 
to the Ising model (referred to as free spins) . The second 
part W num of the coarse-graining transformation replaces 
the free spins inside a block with a single, effective free 
spin. The resulting effective spin model contains con- 
strained spins interspersed with free spins. 
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and ground state fidelities, see Sect. VI 



FIG. 2. Accuracy in the ground state energy of the deformed 
toric code in a 4 x 4 lattice, as a function of the refinement 
parameter \ ( see Sect. VI I, for magnetic field h x — 0.1 in H^ c 
of Eq. |11| On the 4x4 lattice we obtain the exact ground 
state energy by directly diagonalizing the full Hamiltonian 
and use it to asses the precision of the results as a function 
of the computational cost. The computational cost grows 
monotonically with \. Both curves correspond to using a tree 
tensor network (see Sect. W\ as a means to numerically coarse- 
grain the system. However, the lower curve was obtained by 
first applying Wcxact (in this work we refer as the hybrid ansatz 
to the combination of W cxa ct and a tree tensor network, see 
Sect. \Vf . For a fixed value of Xi which roughly corresponds to 
the same computational cost, the hybrid ansatz leads to about 
four more digits of accuracy in the ground state energy. 



Importantly, the local Z2 symmetry, which acts on the 
constrained spins only, is exactly preserved during the 
coarse-graining. In addition, thanks to applying duality 
mappings to the Ising model only within each block of 
spins, the coarse-graining transformation is completely 
local: any local operator of the original model is mapped 
into a local operator of the effective model. Finally, com- 
position of coarse-graining transformations reduces a fi- 
nite system to a small number of spins, which can be 
addressed with exact diagonalization. As it is costumary 
in tensor network algorithms, a variational ansatz for e.g. 
the ground state of the model is then obtained by regard- 
ing the coefficients that characterize the coarse-graining 
transformation as variational parameters. 

Thus, one of the highlights of the approach is that 
it produces a tensor network representation of the wave 
function of the ground state (and low energy states) of 
the model, from which the expectation value of arbitrary 
local observables can be computed. We notice here that, 
depending on the choice of transformation W num , namely 
depending on whether W aum incorporates disentanglers 
or not, our approach can consider arbitrarily large sys- 
tems or is restricted to small lattices. In the first case, one 
can study the renormalization group flow to a fixed point 
(in preparation). In the second case, on which we will 
concentrate here in addition to local observables one can 



As it is costumary in most tensor network algorithms, 
entanglement is central to the present discussion. The 
more entangled a system is, the more costly it is to simu- 
late it with a tensor network ansatz. An important aspect 
of our work is that it highlights the potential role of dual- 
ity transformations in transforming a strongly entangled 
ground state into a weakly entangled one, which is then 
suitable for tensor network algorithms. In particular, the 
toric code for small magnetic field has a robustly entan- 
gled ground state that is mapped into an Ising model with 
large magnetic field, whose ground state is only weakly 
entangled, see Fig. [I] A key of our approach is then to 
be able to map the model to its less entangled dual while 
preserving locality. This is precisely the role of Wcxact- 
After the entanglement in the ground state has been re- 
duced significantly, as illustrated in Fig. [I] W num can 
proceed to coarse-grain the system by discarding high 
energy degrees of freedom. The reduction of entangle- 
ment obtained with Wcxact significantly decreases com- 
putational costs. As a result, with a fixed computational 
cost (as parameterized by some refinement parameter \, 
see Sect. |VI|), use of W cxact produces an improvement of 



four orders of magnitude in the estimate of the ground 
state energy in a 4 x 4 lattice, as shown in Fig [2] As a 
side remark the entanglement structure of the same lat- 
tice gauge theory we consider here has also been studied 
in Ref. [T7D1 



II. THE DEFORMED TORIC CODE 

In this section we briefly review the toric code model, 
an exactly solvable model proposed by Kitaev in Ref. 
11561 We also review a deformation of the toric code 
model obtained by adding a magnetic field on the x di- 
rection, as analysed by Trebst et al. in Ref. 11611 and by 
Hamma et al. in Ref. 11621 for which an exact solution no 
longer exists. The deformed toric code model has a local 
Z 2 symmetry and will be used throughout this work to 
illustrate the proposed coarse-graining transformation. 



A. Toric code 

We consider the toric code modeP^ on a square lat- 
tice £ made of L x L sites and with periodic boundary 
conditions. Recall that in this model a spin-1/2 degree of 
freedom sits on each of the 2L 2 links of the lattice, with 
total vector space 



'tc = 



n«)2L 2 



and Hamiltonian 



A" 



(1) 



(2) 



with J e , J m > 0. Here the star operator A s — \\jes a j 
acts on the spins adjacent to site s and the plaquette 
operator B p = Yl je <Jj acts on all the spins surrounding 
plaquette p, where o~ x and <j z are Pauli matrices, see Fig. 

m 



1. Ground states and topological sectors 

All operators A s and B p commute with each other and 
the ground state subspace of Hamiltonian -f/rc consists 
of the states |£) € Vtc that simultaneously fulfill the star 
constraints 



as well as the plaquette constraints 

B p \0 = \0, Vpe£, 



(3) 



(4) 



and thus have energy —L 2 (J e + J m ). Notice that, on the 
torus, there are L 2 star operators A s and L 2 plaquette 
operators B p fulfilling 



n^- 



sec 



n b p 

pec 



(5) 



Therefore Eqs. [3]and[4]only represent L 2 — 1 independent 
constraints each, that is, a total of 2L 2 — 2 constraints 
in the space of 2L 2 spin-1/2 sites. As a result, there 
are 2 2 = 4 linearly independent ground states. Let us 
introduce the operators 



*i - n "i x *=u 



(6) 



j'eci 



j£c 2 



where C\ and C2 denote the two non-contractible cuts of 
Fig. [3] Operators X\ and X 2 commute with each other 
and with .Hrd and cannot be expressed as products of 
A s '$ and B p 's. The eigenvectors \<b VltV2 ) of Xi and X 2 , 

X l\®v u v 2 ) = Vl\$ VuV2 ), X 2 \$v u v 2 ) = V2\$v u v2), (7) 

where V\,V2 = ±1, form a basis of the ground state sub- 
space. We refer to {|$ +i+ ), | $+,_), !*-,+>, I*-,-)} as 
the ground states of the four topological sectors of the 
model, and we label each topological sector by the pair 
(vi,v 2 ). 

For later reference, we also introduce two operators 



z^U 



z. 



jeh 



n 



(8) 



where l\ and l 2 denote the two non-contractible loops of 
Fig. [3} 
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FIG. 3. In the toric code model, spin 1/2 degrees of freedom 
sit at the edges of a square lattice C. Here we consider a lat- 
tice C, with periodic boundary conditions on both directions 
(torus) . A star operator A 3 acts on the four spins surrounding 
site s £ C, whereas a plaquette operator B p acts on the four 
spins spins surrounding a plaquette p. The non-contractible 
cuts ci and ci and the non-contractible loops l\ and h are the 
support of operators X\ and X2 in Eq. [6] and of operators Z\ 
and Z2 in Eq. fsj respectively. 



electric charge on site s if A s |£) = — 1£). Similarly, we 
say that state |£) contains a magnetic vortex (or magnetic 
monopole) in plaquette p if B p \£) — — 1£). 

Operator erj acting on a ground state l*^,^) produces 
an excited state with a pair of electric charges sitting on 
the two sites s and r connected by link j. Indeed, since 
for site s (equivalently, for site r) we have A s o~ z : 
it follows that 



-a]A s 



A s (a 



j |^«l,«2 



)) 



-o-*-A s \$ VuV2 ) 



[<rj\$v u v 2 



))■ (9) 



State <Tj\& Vu v!t) has an energy 4J e above the ground state 
energy, since each violation of a star constraint increases 
the energy by 2J e . More generally, any other state that 
fulfils Eqs. [3] and [4] except in two sites, again correspond- 
ing to the presence of two electric charges, has energy 
4J e above the ground state. This implies that a pair of 
electric charges created locally by acting with er| on the 
ground state can be separated in space without a change 
in energy. We say that electric charges are deconfined. 
On the other hand, operator o~ x acting on a ground state 
\^ Vl ,v 2 ) produces a pair of magnetic vortices (monopoles) 
sitting on the two plaquettes p and q that contain link 
j. For plaquette p (equivalently, for plaquette q) we have 
BpG? = —o-^Bp, and therefore 



2. Electric charges and magnetic vortices 

Excited states of the toric code Hamiltonian i?rc are 
characterized by violations of the star and plaquette con- 
straints, Eqs. [3] and El We say that state |£) contains an 



B,, 



o- x \$ 



))=-o-*B p \$ VuV2 )=-(a*\$ Vl , V2 )). 

(10) 

1 has an energy 4J m above the ground 
state energy. Any other state with just two magnetic 
vortices has energy 4J m and we say that magnetic vor- 
tices are also deconfined. 



State <r*|$ vl , Ua 



B. Hamiltonian deformation of the toric code 

The toric code Hamiltonian ifrc m Eq. [2] is exactly 
solvable. Here we will be interested in a non-solvable 
deformation of the toric code ob tained by introducing a 
magnetic field in the x direction] 161 * 162 ! 

H* c = -j e j2A s -j m J2 B p- h *I2 a J- ( n ) 



The magnetic field lifts the ground state degeneracy. No- 
tice, however, that since y\. a? also commutes with the 
operators X\ and X2, Hamiltonian -ff^c s ^ m decomposes 
into four sectors (vi, ^2). The ground state of any of the 
sectors may no longer fulfill the plaquette constraints, 
indicating the presence of pairs of magnetic vortices, for 
which J^ . o X j also acts as a hopping term. 

The limit J e 3> J m ,h x , where low energy states fulfill 
the star constraints of Eq. [3] an d th erefore contain no 
electric charges, is well un derst oocP^ (it corresponds to 
the Z2 lattice gauge theoryLUiiJ which in turn is dual to 
the quantum Ising model, as reviewed in Sect. Ill I. In 



this regime, the ground state phase diagram of H^ c con- 
tains two phases, one in which pairs of magnetic vortices 
are deconfined, and another with a Bose condensate of 
magnetic vortices. 



1. Deconfined phase 

For small values of the magnetic field h x , h x <C 1, 
the model is in a deconfined phase, in which a pair of 
(dressed) magnetic vortices created locally can be sep- 
arated an arbitrary distance incurring only a finite en- 
ergy penalty. [Notice that in this work the term 'decon- 
fined' refers to magnetic vortices, and not to the electric 
charges, which by construction are not present in the low 
energy sector of the model.] 

The deconfined phase is a topologically ordered phase 
with four nearly degenerate ground states on the torus, 
one for each sector (ui,i>2). The sector (+, +) has the 
smallest energy and the energy separation A to another 
sector vanis hes exponentially fast with the linear size L 
of the lattice^, 



A 



,-L/i 



(12) 



where £ is a finite length scale that vanishes for h x = 0, 
i.e. A = for the undeformed toric code. 

To further c haracterize this phase, let us introduce the 
string operatoJ^ 54 1 1711 ^^ 



x^U 



(13) 



j£c 3 



where C3 is the cut of Fig. H connecting a pair (po,Pi) 
of plaquettes that are as distant as possible in the torus. 
Since X 3 anticommutes with B po and B pi , and commutes 
with the rest of plaquette terms, an argument similar to 



that in Eq. 10 above shows that this operator acting 
on a ground state \& Vl ,v 2 } f° r h x = (undeformed toric 
code) produces a state with a pair of magnetic vortices / 
monopoles sitting on plaquettes po and p\ . This state is 
orthogonal to the ground state and therefore, for h x = 0, 
the ground state expectation value of X3 vanishes, 



(*3> = 0. 



(14) 



In the thermodynamic (or large L) limit, (X 3 ) vanishes 
for the whole deconfined phase and it is used as a (non- 
local) disorder parameter. 

Alternatively, the deconfined phase can also be charac- 
teriz ed b y the scaling of the expectation value of Wilson 
loopsSS. For every contractible loop I on C, see Fig. |4j 
we can define a Wilson loop operator Z[l] by 



z[i] S n°?- 



(15) 



Operator Z[l] amounts to the product of operators B p 
for all plaquettes p contained in the interior of loop I. 
Therefore, for any ground state |$ Wl ,t) 2 ) for h x = (un- 
deformed toric code), we have the expectation value 



(Z\. 



1. 



(16) 



More generally, in the deconfined phase Wilson loops 
obey a perimeter lav)^* , in the sense that they decay as 



-ap(l) 



(17) 



where p(l) is the length of the loop I, and where a > 
vanishes for h x = 0, recovering Eq. [16] for the unde- 
formed toric code. The behaviour of the Wilson loop 
with respect to its length is related to the dependence 
of the potential between two static electric charges with 
their mutual distance, which is in turn induced by their 
interaction with the gauge fields. In the deconfined phase 
this potential is such that the force between the two 
charges does not increase with their distance. This trans- 
late into a decay of the Wilson loop proportional to the 
perimeter of the loop, as can be f ound in standard text- 
books about lattice g auge theor iea^^ and in the original 
papers on the subjedP^EHEZI 



2. Spin polarized phase 

For large values of h x , the model is in a state where 
the spins are polarized in the x direction. 

This phase has a unique global ground state, corre- 
sponding to the ground state of sector (+,+), and the 
energy separation to the other sectors grows linearly in 
the system size L, 



L. 



(18) 



The spin polarized phase can be interpreted as a Bose 
condensate of magnetic vortices /monopoles, and is char- 
acterized by a non-vanishing [actually, positive] value of 
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FIG. 4. Cut C3 connects two plaquettes po and pi that are 
as far apart in C as possible. Example of loop / that has a 
length of eight spins, |/| — 8, and encloses an areas of four 
plaquettes, E(7) — 4. 



the (non-local) order parameter (X3), 

(*3, * 0, 



(19) 



which can be interpreted as the square of the expectation 
value of a creation operator for a single magnetic vortex 
(see Sect. HID), with (X3) = 1 for the completely polar- 
ized state at h x — ¥ 00. Alternatively, it is characterized 
by an area law of the ground state expectation value of 
Wilson loops, 



,-/9o(0 



(20) 



where a(l) is the area enclosed by the loop I, and where 

p — 00 for h x — > 00, in which case (Z[t\) = 0. 



3. Continuous quantum phase transition 

For an intermediate value of h x of about h x = 
0.3283221, the system undergoes a continuous quantum 
phase transition. This is in the same universality class 



of that of the 2D quantum Ising model (see Sect. HID) 



or the 3D classical Ising model, but it is special in that 
it separates two phases which cannot be distinguished 
by a local order parameter. Instead, the transition can 
be characterized by the non-local order parameter (X3), 
which detects the formation of a Bose condensate of mag- 
netic vortices /monopoles, or by the scaling of Wilson 
loops as a function of their size, and by the presence / 
absence of an approximate ground state degeneracy. 



Apart from the Z 2 x Z 2 symmetry generated by oper- 
ators X± and X2, H^ c is also invariant under unitary 
conjugation by any star operator A s , 



A S H TC A S — H TC , 



Va € C, 



(21) 



which follows from recalling that [H^ c , A s ] — and from 
noticing that A s = A\ = Aj 1 . Therefore, the deformed 
toric code Hamiltonian H^ c has a local Z2 symmetry, 
generated by unitary transformations of the form 



A s 



'TC 



'TC> 



(A s y 



1. 



(22) 



which simultaneously flip the four spins included in a star 
s, for any choice of star s £ C. 

This local symmetry implies the presence of L 2 — 1 
constants of motion, namely the eigenvalues ±1 of L 2 — 1 
independent star operators A s . In particular, a state that 
fulfils all the star constraints of Eq. [3] A s \£) = |£), that 
is, a state that is invariant under the local symmetry, 
remains so under a time evolution generated by H^ c . It 
is important to note, however, that the vector space Vtc 
also contains states that are not invariant under A s . For 
instance, we have already seen that a state with a pair 
of electric charges, Eq. [9j violates the star constraint in 
two sites, where it transforms as A s \£) = — 1£). 



D. Entanglement 

Entanglement is another aspect of the deformed toric 
code model that will play a key role in for subsequent 
discussions on how to coarse-grain the system. 

When the lattice is in a pure state |\1/) £ Yc, the en- 
tanglement between a block of spins and the rest of the 
system can be measured by means of the Von Neumann 
entropy S of the reduced density matrix p of the block of 
spins, known as entanglement entropy, 



S = -tr (p\og 2 p) , p = tw|WI 



(23) 



Fig. [T] shows the entanglement entropy of a block of 
four spins as a function of the magnetic field h x , when 
the system is in the ground state 1$+,+) of H^ c . No- 
tice that for h x = 0, corresponding to the (undeformed) 
toric code, the entanglement entropy of the block is sig- 
nificantly large, indicating that the ground state is very 
entangled. This is in sharp contrast with the entangle- 
ment entropy for large h x , which tends to zero. In this 
second case, each spin is completely polarized in the x 
direction and the ground state is a product, unentangled 
state. 



C. Local symmetry 



III. Z 2 LATTICE GAUGE THEORY AND THE 
QUANTUM ISING MODEL 



The symmetries of the deformed toric code Hamilto- 



nian H^ c in Eq. 



In this section we briefly review Z 2 lattice gauge 
11 play a central role in this paper. theory^ 5 ! 152 ! 154 * 1 -^}!!!^) j n two spatial dimensions and 



its c onnection to the two-dimensional quantum Ising 
model^ 4 * 181 * 182 ! . We also introduce a convenient graph- 
ical notation for the different Hamiltonian terms in the 
Hamiltonian, and pieces of the stabilizer formalism, in- 
cluding the transformation rules of Pauli matrices un- 
der so-called CNOT gates. The stabilizer formalism pro- 
vides us with the natural language to describe the coarse- 
graining transformation of systems with a local symme- 
try, as described in Sect. [V] 



A. Zi lattice gauge theory 

The considered restriction J e 3> J m ,h x on H^ c en- 
sures that low energy states fulfill the star constraints 
of Eq. [3J If, instead, we impose those constraints by 
truncating the vector space of the model, we are left 
with the vector space Vlgt C V tc of Z 2 lattice gauge 
theor ^ 115 ' 152 ' 1 ^ , 

Vlgt s{ |0 e Vtc : ^.10 = 10 Vs e £} (24) 

i>L 2 + l 



= (C 2 f 



(25) 



When projected onto Vlgt, Hamiltonian -ff^c m Eq- HI] 
becomes 



-ffLGT — — 2_^i Bp ~ hi 



2^ U 3 



(26) 



where we have neglected a constant term —L 2 J e and we 
have set for simplicity, without loss of generality, J m = 1. 
It is also understood that each term in i?LGT is pro- 
jected on to Vlgt- This model has been extensively stud- 
ied before^- 2 * 154 * 179 * 1 ^. It was originally introduced as a 
model that presents two different phases that cannot be 
distinguished with a local order parameter. By increasing 
the intensity of the magnetic field one can indeed drive a 
transition from the deconfined phase (at small magnetic 
fields) to a spin polarized phase (at large magnetic field). 
The phase transition around h x ~ 0.3 is in the 3D Ising 
model universality class. The Z 2 lattice gauge theory 
has also been intensively studied because it presents the 
ideal testing ground for new ideas, since its relative sim- 
plicity (as compared to gauge theories with larger gauge 
groups) allows to obtain v ery precise numerical results 
in Monte Carlo simulations^ 8 -^. Further interest in the 
model comes from a conjecture that relates its critical 
properties with the ones of the finite temperature decon- 
fining tran sitio n of an SU{2) lattice gauge theory in 3 + 1 
dimensional 8 -^. Finally, the existence of a duality trans- 
formation between the Z 2 lattice gauge theory and the 
quantum 2D Ising model, as we review in Sect. HID 



has inspired a lot of work in finding the relevant order 
parameter for confinements'. 

By construction, the Z 2 lattice gauge can be regarded 
as a low energy, effective model of H^ c (for J e 3> J m , h x ) 
and therefore has the same ground state phase diagram. 

Hamiltonian i?LGT has a Z 2 x Z 2 symmetry corre- 
sponding to operators X\ and X 2 (properly projected 
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FIG. 5. (i) Graphical representation of the Hamiltonian .Htc 
for a 3 x 3 lattice C. A four-spin star operator A a acts on each 
site s £ £., and a four-spin plaquette operator B p acts on each 
plaquette p £ C. (ii) Hamiltonian Htc can also be defined on 
more general lattices. The example shows an irregular lattice 
with a five-spin star operator A a and a three spin plaquette 
operator B p . (iii) Additionally, we represent the magnetic 
field of the deformed toric code Hamiltonian H TC by means 
of a yellow shade surrounding each qubit on which —h x a x 
acts. 



onto Vlgt)- However, it is worth emphasizing that, in 
contrast with Hamiltonian -H"xc> Hamiltonian H^ct does 
not have a local symmetry. Indeed, H^ct is defined on 
the subspacc Vlgt <= Vtc oi Eq. [24j which is made 
of vectors that fulfill the star constraints. That is, each 
star operator A s acts as the identity operator in Vlgt, 
and therefore the assertion that i?LGT is invariant under 
transformation A s is an empty statement. 

The present work aims at developing a coarse-graining 
scheme for lattice models with a local symmetry. For con- 
creteness, we will most of the time restrict our attention 
to the subspace of states that fulfill the symmetry con- 
straints, that is, Vlgt- In this case, the coarse-graining 
transformation can be used also to obtain a low energy, 
effective description of lattice gauge models, as illustrated 
here in the context of Z 2 lattice gauge theory. 



B. Graphical representation of Hamiltonians 



In this paper we will represent star and plaquette oper- 
ators graphically by just drawing the corresponding lat- 
tice. For instance, Fig. [5h) represents the toric code 
Hamiltonian Htc of Eq. [2J The toric code Hamiltonian 
Htc can be defined on a more general lattice by consid- 
ering star operators A s and plaquette operators B p that 
involve a variable number of spins or qubits, see Fig. 
[5jii). In that case, we will also often just denote Htc by 
drawing the underlying lattice. 

In addition, the presence of a magnetic field or, more 
generally, of related additional interactions, will be rep- 
resented by a yellow shaded region enclosing all relevant 
qubits. For instance, Fig. pniii) represents the deformed 
toric code Hamiltonian i?xc °f ^q- 11 



C. Qubits and CNOTs 

In order to describe our coarse-graining scheme for the 
deformed toric code model and the Zi lattice gauge the- 
ory, it is convenient to introduce first some basic pieces of 
nomenclature and formalism frequently used in the area 
of quantum information. 

A spin- 1/2 degree of freedom is an example of a quan- 
tum bit or qubit. We will work with two preferred basis 
of a qubit. One basis corresponds to the eigenvectors |0) 
and |1) of <r z ; the second one, to the eigenvectors |+) and 
\-)oia x , 



<7*|0) = |0), 



o*|+) = |+), (27) 

°*\-) = -h), (28) 



where 



1+) = ^, H = ^^. (29) 



V2 ■ 



V2 



A control-not gate or CNOT gate is a unitary trans- 
formation Uqnot on two qubits given by 

t/ CNOT = |0)(0|®I + \l)(l\®a x (30) 

= ISl+X+l + cr- O J— X— U (31) 

where the first and second qubits are referred to as the 
control and targ et q ubits, respectively. We use the sta- 
bilizer formalism!^ to study how operators evolve under 
the application of a CNOT gate. Specifically, under con- 
jugation by a CNOT, Pauli matrices on the control and 
target qubits are transformed according to 




I <g> a z <H> a z ®a z , 


a z <g> I «-> a z <g> I, 


(32) 


1(g) a x ^ l®a x , 


a x <E) I «-» a x <E) cr x . 


(33) 



FIG. 6. (i) A CNOT is represented with an arrow from the 
control qubit to the target qubit. When applied on two ad- 
jacent qubits, it reconnects the target qubit. This affects the 
two stars r, s to which the control qubits belongs, which be- 
come r' and s' , as well as the two plaquettes p, q to which the 
target qubit belongs, which become p' and q' . (ii) A sequence 
of CNOTs can also be used to focus a star operator A s on 
a single qubit, on which it acts as A s i — crj. As a result, 
the neighbouring site r increases the number of qubits con- 
nected to it, becoming r' . Notice that the single-qubit star 
can be represented inside any of the plaquettes surrounding 
site r' . In this work we restrict our attention to the subspace 
of states that are invariant under star constraints, Eq. [3] The 
single-qubit star A s i forces the qubit to be in state |+), and 
therefore unentangled with respect to the rest of the qubits. 
We will use this fact to remove it from the effective description 
of the system at low energies, (iii) A sequence of CNOTs can 
be used to focus a plaquette operator B p on a single qubit, 
on which it acts as B p i = erf/. As a result, the neighbouring 
plaquette q expands to q' . Notice that the single-qubit pla- 
quette can be represented as connected to any of the sites of 
plaquette q' . 



Given the toric code Hamiltonian on a given lattice, 
the action of a CNOT between two neighbouring qubits 
is to reconnect the target qubit according to Fig. [6Vi) . 
Notice that we represent a CNOT gate by an arrow from 
the control qubit to the target qubit. 

Following Ref.1££l, a sequence of CNOTs involving the 
qubits of a star can be used to reduce the star operator 
to a single qubit, Fig. p[ii). The star operator reduces to 
a x acting on that qubit. Similarly, a sequence of CNOTs 
involving the qubits of a plaquette can then be used to 
reduce the plaquette to a single qubit, as illustrated in 
Fig. [6]jiii). The plaquette operator reduces to a z acting 
on that qubit. 

Of course a CNOT between two qubits will not only 
modify the star and plaquette operators. For instance, 
according to Eq. [33j a magnetic field a x acting on the 
control qubit will be transformed into a a x <g a x interac- 
tion between control and target qubits, see Fig. [7j 

In the rest of the paper we refer to qubits as spins. 



D. Duality transformation: quantum Ising model 

The Z 2 lattice gauge theory has plaquette degrees of 
freedom that can take two values, namely the eigenvalues 
±1 of the plaquette operators B p , corresponding to the 
absence/presence of a magnetic vortex on that plaquette. 
Thus plaquettes behave as spin-1/2 degree of freedom. In 
addition, as argued in Eq. 10 the operator a x can flip 



the value of the two nearest neighbour plaquettes p and 
q that share the link j. This interpretation of the Z 2 
lattice gauge theory as a spin model can be materialized 
explicitly by means of a well-known duality transforma- 
tion that maps it into t he q uantum Ising model with a 
transverse magnetic fielcP^l ( se e appendix \M , 
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FIG. 7. The magnetic field a x acting on control qubits ex- 
pands into a multi-qubit interaction a x ® o x d§ ■ ■ ■ including 
all target qubits. (i) In concentrating a star operator on a 
single qubit, the magnetic field on that qubit spreads to all 
the qubits originally involved on that star. In the subspace 
of locally symmetric states, Eq. [3l we can remove the con- 
trol qubit (now in state |+) and thus unentangled from the 
rest of the system, see caption of Fig[6[ii)) from the effective 
description, (ii) In concentrating a plaquette operator on a 
single qubit, the magnetic field a x of the rest of qubits on 
the original plaquette turn into a a x <g> o x interaction between 
those qubits and the target qubit. 
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FIG. 8. (i) Left: graphical representation of the toric code 
Hamiltonian Htc, equivalently, -Hlgt for h x — 0, for a lattice 
C made of 4 x 4 sites or 4 2 x 2 = 16 spins. Right: graphical 
representation of the dual quantum Ising Hamiltonian flising- 
There are 4x4 — 1 = 15 spins with a single-spin plaque- 
tte operator, that is, with a transverse magnetic field — J m /x 2 . 
Notice the vacancy and two additional topological spins which 
are not subject to any star or plaquette constraints. Opera- 
tors X\ and X2 act each on one of these qubits. 



where the spins are placed on the sites of a L x L lattice 
£duai dual t t ne original lattice C (i.e., they can be iden- 
tified with the plaquettes of C) , and fi x and [i z are Pauli 
matrices, see Figs. [8][9] However, the resulting Ising 
model inherits two unconventional elements. On the one 
hand, lattice £ dual has one vacancy, i.e. there are only 
L 2 — 1 spins, reflecting the fact that the Z2 lattice gauge 
theory only had L 2 — 1 independent plaquettes. [It is 
worth noticing that the /z* terms around the vacancy ex- 
plicitly break the Zi symmetry of the Ising model.] On 
the other hand, there are two additional non-local spin- 
1/2 degrees of freedom, associated to operators X x and 
X2 , each of which is coupled to L spins along one of the 
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FIG. 9. (i) Graphical representation of the Hamiltonian H^ c , 
equivalently Hlgt, including the magnetic field —h x ^2,o x . 
(ii) The duality transformation maps the magnetic field to a 
nearest neighbour ferromagnetic Ising interaction — h x fi x ®u x . 
For simplicity, transverse magnetic field terms —J\jl z of -Hi s j ng 
are not depicted here. Notice that the four spins surrounding 
the vacancy, denoted U, D, L and R, as well as all boundary 
spins, follow a different interaction pattern. On the one hand, 
U and L have one interaction term less than the rest of the 
spins, but have instead an additional term —h x ^i x acting on 
them, (iii) On the other hand, the nearest neighbour inter- 
action between spins p and q at both sides of a boundary is 
mediated by a coupling ~h x fip (8) X a ® fi x that includes one 
topological spin (here, X a corresponds to one of the non-local 
operators X\ or X2 defined in Eq. pi. In the sector (+,+), 
the topological spins are in state |+)<8)|+) and boundary spins 
are coupled with a regular Ising interaction —h x ^L x ®^i x , corre- 
sponding to periodic boundary conditions (BC) in both direc- 
tions. However, in e.g. sector (+, — ) the topological spins are 
in state |+) ® |— ) and the coupling between the spins in the 
uppermost and lowermost rows becomes antiferromagnetic, 
h x fj,p ® /Xq, corresponding to antiperiodic BC in the vertical 
direction, while in the horizontal direction the BC are still 
periodic. In general, sector (vi,V2) the BC in the horizontal 
(respectively vertical) direction will be periodic/antiperiodic 
depending on whether V\ (respectively V2) is +/— . 



two directions in £ dual , see Fig. [9J 

By ignoring the presence of the vacancy and non-local 
spins, which do not affect local observables in the ther- 
modynamic limit, the ground state phase diagram of 
Hlgt (and thus of H^ c for J e 3> J m , h x ) can be recov- 
ered fro m the well-know n ground state phase diagram 
of His . n mMMMMEMM by inver ting the duality trans- 
formation. Recall that the 2D quantum Ising model at 
zero temperature has a disordered phase for small h x 
and an ordered phase for large h xi as characterized by 
the vanishing (respectively non-vanishing) value of the 



10 



spontaneous magnetization m. 

1 



(i) 



■in. 



EW. 



(35) 



which is a local order parameter, and that the two phases 
are separated by a continuous q uantum phase tra nsition 
that occurs, indeed, at h x « o.3 3 m ' 144 l lii7 m m™ 

In particular, we can gain insight into the meaning of 
the nonlocal order parameter (X 3 ) of the deformed toric 
code by noticing that operator A3 in Eq. [13] can be 
written as 



X 



3 — H'poH'p!) 



(36) 



where po and p\ are two spins in £ dual (or plaquettes 
in C) separated by O(L) sites. In the limit L ^S> 1 of a 
large lattice, where {n Po fi> Pl ) is expected to factorize into 
(jUp )(/Up 1 ) since plaquettes po and pi are far apart, this 
expectation value equals the square of the spontaneous 
magnetization m x , 

(X 3 ) = (n x P0 H x Pl )^(n x P0 )(n x pl ) = m 2 x (37) 

Thus, we can interpret (A3) as the square of the expecta- 
tion value of a (fictitious) creation operator of single mag- 
netic vortices, and interpret a vanishing/non-vanishing 
value for (X3) as indicating the absence/presence of a 
Bose condensate of magnetic vortices. 



IV. ENTANGLEMENT RENORMALIZATION 

In this section we briefly review one possible route to 
coarse-graining a lattice mo del in real spa ce based on en- 
tanglement renormalizatiort 1 - 41 * 142 * 145 * 1 ^ 1 -*. Entanglement 
renormalization is a speci fic re al-space implementation of 
the renormalization group^l. The coarse-graining trans- 
formation can be used to de fine t he multi-scale entangle- 
ment renormalization ansatz^^ (MERA), a variational 
ansatz for ground states and low energy states of a lo- 
cal Hamiltonian. We also review how to use entangle- 
ment renormalization and the MERA in the presence 
of a global internal symmetry, in such a way that the 
symmetry is exactly preserved and exploited to reduce 
computational costs. Then, in Sect. [Vl we will explore 
how to generalize entanglement renormalization and the 
MERA in the presence of a local symmetry. 



A. Coarse-graining transformation 

Given a lattice model, characterized by a lattice C and 
a Hamiltonian H , a coarse-graining transformation aims 
to produce a new, simplified lattice model, characterized 
by a coarse-grained lattice £ with fewer sites and an 
effective Hamiltonian H' . Here we are interested in a 
coarse-graining transformation W 
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FIG. 10. (i) The coarse-graining transformation W defines 
a linear map between the space Y& of lattice C' and the 
space Vc of lattice C. (ii) Isometric transformation W that 
decomposes as a product of disentanglers u and isometries w. 
(iii) Constraints fulfilled by disentanglers and isometries, see 
Eq. HO] 
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FIG. 11. Coarse-graining of a local operator o acting on three 
contiguous sites of C, into a local operator o' acting on three 
contiguous sites of C' . (i) Diagrammatic expression for o = 
W^oW, where W decomposes as a product of disentanglers 
u and isometries w. (ii) Using that u'u = /, where / is the 
identity on two sites of £, most disentanglers can be removed, 
(iii) Using that w^w — I, where / is the identity on one site 
of £', most isometries can be removed, so that (iv) o' acts as 
the identity in all but three sites of C ' . 



that defines a linear map from the space Vc of the orig- 
inal model to the space Vc of the effective model, see 
[Notice that, for consistency with Ref. I142[ 

to V c . } 



Fig. 10 ;i) 



we actually regard IT as a map from Vc 

More specifically, we consider a coarse-graining trans- 
formation W that maps local operators in L (that is op- 
erators that act non-trivially on a finite subset of neigh- 
bouring sites of £) into local operators in C! . This prop- 
erty is automatically fulfilled if W is an isometric tensor, 



W^ : V c -> V c > 



(38) 



W^W 



(39) 
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FIG. 12. Effective Hamiltonian H' = W 1 HW expressed in 
terms of the original Hamiltonian H and disentanglers u ans 
isometries w. If H is a sum of local terms, then H' can also 
be expressed as a sum of local terms, see Fig. [TTJ 



that decomposes as the product of isometric tensors u 
-(in), 



and w, see Figs. 10 
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V/U = /, VTW = I 
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(40) 



known as disentanglers and isometries, which corre- 
sponds to one step of entanglement renormalization^^. 
Indeed, as shown in Fig. [TT] in this case local operators 
are mapped into local operators. In particular, if the 
Hamiltonian H of the original lattice model decomposes 
as the sum of local terms, then the effective Hamiltonian 
H' , depicted in Fig. 12 will also decompose as a sum of 



local terms. 

The examples of Figs. [10][T6| used to illustrate this 
section correspond to one dimensional lattices for sim- 
plicity. Analogous constructions for two dimensional lat- 
tices span three dimensions, see e.g. in Ref. 11441 and 
11451 Their graphical representation is significantly more 
involved and is not required in order to introduce the 
basic elements necessary for the present discussion. We 
postpone until Figs. [20] and [23] of Sect. [V] the explicit 
representation of three dimensional structures. 



B. Renormalization Group flow 

Successive applications of coarse-graining transfor- 
mations {W, W, W", ■ ■ ■ } produce a sequence of in- 
creasingly coarse-grained lattices {£,£',£",■■■} to- 
gether with a sequence of effective, local Hamiltonians 
{H, H' ,H" , • • • }, see Fig. 13 If at each step the coarse- 



graining transformation projects onto the low energy sub- 
space of the Hamiltonian, then the sequence of effective 
Hamiltonians defines a discrete renormalization group 
(RG) flow towards a fixed point model that captures the 
low energy/large scale properties of the original model. 

Suppose now that the original lattice £ is made of a 
finite number N of sites. Then the effective lattice £' 
is made of a smaller number N' = N/b of sites, where 
b = 2 in the example of Fig. [lO^ii). In particular, after 
0(log b (N)) applications of the coarse-graining transfor- 
mation, lattice £ is reduced to a small lattice £ top whose 
size is independent of N, and such that H top is amenable 
to exact diagonalization techniques. Notice that a state 




£" 
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FIG. 13. By composition, coarse-graining transformations 
{W, W' , W", • • • } produce a sequence of increasingly coarse- 
grained lattices {£,£',£.",■■■} with effective Hamiltonians 
{H,H',H",.-.}. 



|\E» top ) £ V£to P depends on a small number of parameters 
that is also independent of the original system size N. 



C. MERA 

A state |<I/ top ) £ V^top of this reduced lattice can be 
mapped into a state | VP) £ Vc of the original lattice by re- 
versing the coarse-graining transformations, see Fig. [TT] 



|*) = WW'W" ■■■\^ top ). 



(41) 



This means that the top vector |^ top ) together with the 
sequence of transformations {W, W' , W" , ■••}, as char- 
acterized by a sequence of disentanglers and isometries 
{{u, iv}, {u',w'},{u", w"}, ■ ■ ■ }, can be used as an effi- 
cient representation a state \^f) £ Vc- This represen- 
tation is the multi-scale entanglement renormalization 
ansatz (MERA) , and can be used as a variational ansatz 
for e.g. the ground state of H. This is achieved by opti- 
mizing the coefficients in {{u, w}, {u 1 , w'}, {u", w"}, ■ ■ ■} 
and |v[/ top ) so as to minimize the expectation value 
( 1 $ r |.ff|\I< f ) using the techniques discussed in Ref. 11451 
More generally, by replacing |vI/ to P) with an isometry 



W top that defines a x* -dimensional space 



'X > 



a whole 



subspace V* C Yc can be represented. If \a) denotes an 



orthonormal basis in 



then 



\^ a ) =WW'W"---W top \a), a = l,---,x* (42) 

is an orthonormal basis in V* C V^. In this way, the 
MERA can be used to approximate e.g. a low energy 
subspace, including the ground state(s) and several low 
energy excited states of H, as we will do in Sect. |VI| 

When the state |\I/) in Eq. fid] is translation invariant, 
which is the case we have discussed so far, the MERA 
is specified with 0{logb{N)) parameters, namely those 
that characterize the 0(log b (7V)) pairs of disentangles 
and isometries {{u,w},{u' ,w'},{u" ,w"},- ■ ■}. Thus, 
the MERA offers an efficient description of certain states 
of the lattice model, whose vector space V^ has dimen- 
sion 0(exp(JV)). (More generally, a generic MERA is 
specified by N independent disentanglers and isometries, 
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FIG. 14. Example of a state of a lattice C made of 16 sites, 
|*) = WW'W"\y top ), expressed as a MERA made of three 
layers of disentanglers and isometries corresponding to three 
coarse-graining transformations W, W' and W" . 



and therefore depends on O(N) parametera^Sl) . In ad- 
dition, a scale invariant state \^/) G Yc of the form 41 
can be described using O(l) parameters by choosing all 
disentanglers and isometries to be identical. The scale 
invariant MERA is useful to represent the ground state 
of fixed-points of the RG flow , whether corresponding to 
critical s ystemi 1 41 * 142 * 193 -* 1 ^ or to systems with topolog- 
ical ordtnPaSSl 

In Sect. IVII we will use a translation invariant coarse- 
graining transformation for the deformed toric code 
model, which will lead to a translation invariant MERA 
for the ground state(s) and low energy excited states of 
H^ c . However, for simplicity we will not attempt to use 
disentanglers in their full form (see Sect. VF for details) 



as a result of which the computational cost will grow as 
0(exp(ViV)). The limits of small and large magnetic 
field of Harniltonian H^ c , i- e - h x = and h x = oo, 
correspond to fixed points of the RG flow, namely to a 
topologically ordered fixed poinlpSSl anc j to a trivial fixed 
point, respectively. In these two cases, we will be able to 
represent the ground state with a scale invariant MERA. 
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FIG. 15. (i) Invariance of a Harniltonian H under a global 
symmetry implemented by simultaneously acting with one- 
site transformation r on all sites of C, see Eq. 



44 



(ii) Dis- 

entangler u and isometry w invariant under the action of the 
single site transformation r, represented by a circle, acting 
on all upper/lower indices of these tensors. Notice that the 
upper index of the isometry w is transformed according to 
r', which might be a different representation of the symmetry 
group. 
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FIG. 16. Sequence of equalities showing that a global sym- 
metry R = r® N of Harniltonian H is exactly preserved by a 
coarse-graining transformation W that is the product of dis- 
entanglers u and isometries w, provided the latter are also 
invariant under transformation r, Fig. |15| Indeed, the effec- 
tive Harniltonian H' = W^ HW is invariant under the global 



transformation R' — (r 



l\®N' 



D. Global symmetry 

Our goal is to coarse-grain a lattice model with a local 
symmetry. It is instructive to first consider the sim pler 
and better understood case of a global symmetrjEED. 

Let r be a unitary transformation acting on one site of 
lattice C, and let 



R=r 



®N 



(43) 



be the unitary transformation that results from applying 
the same transformation r on all sites of C simultane- 
ously. For instance, in the case where a site is described 
by a spin-1/2 degree of freedom, the one-site unitary 
transformation r could correspond to a Pauli matrix, say 
r = ia x , in which case R — {io x )® N corresponds to si- 
multaneously rotating all spins in the lattice by and angle 
7r in the x direction. 



The lattice Harniltonian H is invariant under R if 

RHR) = H, (44) 

an equation that can be represented diagrammatically as 



in Fig. 15 '}). In this case, we say that the lattice model 
has a global symmetry R. Recall that if R is a global 
symmetry of H, then R 2 and i? _1 are also a global sym- 
metries of H and, more generally, the global symmetries 
of H always form a group. In the case of r = ia x , where 
r 2 = I, this group is Z^. We also notice in passing that 
this particular choice r = ia x is actually a global symme- 
try of the deformed toric code Harniltonian fl^O si nce it 
is included as part of its local symmetry. There, R corre- 
sponds to the product of N/2 star operators A s , chosen 
according to a checkerboard pattern. 

The presence of a global symmetry is a fundamental 
property of a lattice model, since the eigenvectors of H 
are organized according to representations of the symme- 
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FIG. 17. Diagrammatic representation of a Hamiltonian with 
a local symmetry. In this example the local symmetry is im- 
plemented by unitary transformations r <g> r <g> r <g> r , acting on 
some blocks of four contigous sites. 



try group. It is therefore of interest to consider a coarse- 
graining transformation W that preserves the global sym- 
metry of a model. With exception of models with spon- 
taneous symmetry breaking, which will not be considered 
here, the low energy subspace automatically inherits the 
symmetries of a model and, in principle, there should be 
no need to enforce symmetry preservation during coarse- 
graining. However, due to the numerical (and therefore 
approximate) nature of W, the symmetry may be lost 
during coarse-graining, unless some measures are put in 
place to explicitly protect it. 

In the case where W is the product of disentanglers 
u and isometries w, a global symmetry is automatically 
preserved if these tensors are chosen to be symmetric 
themselves, that is 



(r r)u(r ® ry — u, (r ® r)w(r ) 



'M - 



(45) 



where r' is a unitary matrix acting on a site of £' ', which 
in general will be transformed according to a different 
representation of the same symmetry group. The invari- 
ance of u and w implies that the action of the group 
commutes with these tensors, see Fig. Il5jii). As shown 
in Fig. |16[ in this case indeed the coarse-grained Hamil- 
tonian H' = W^HW will also be invariant under the 
global transformation R' = (r')® N . 

Using symmetric tensors u and w does not only ensure 
exact preservation of the global symmetry, but it also 
allows for important computational savings. On the one 
hand, symmetric tensors depend on less parameters than 
generic tensors, leading to a more compact variational 
ansatz. On the other, manipulating symmetric tensors 
(for instance in order to compute the expectation value 
of a local observable) has lower computational cost than 
manipulating generic tensors. 



E. Local symmetry 

Let us now assume that the lattice model has a local 
symmetry. In this case the Hamiltonian H is invariant 
under unitary transformations that act on a small subset 
of neighbouring sites, as exemplified by Eq. 



21 for the 



deformed toric code Hamiltonian H^ c , that is invariant 
under star operators A s acting on four contiguous spins. 
This situation is diagrammatically represented in Fig. 17 



How can a local symmetry be exactly preserved during 
coarse-graining? And how can it be exploited to obtain a 
more compact variational ansatz and to reduce the com- 
putational cost of simulations, as in the case of a global 
symmetry? These questions were already addressed by 
Fradkin et al. in Ref. 11981 and before it by the seminal 
work of Migdal et al. in Ref. 11991 In the next section we 
take a fresh look at them, by applying the ideas of En- 
tanglement Renormalization summarized in Sect. IV to 
this subject. Next section is indeed devoted to describe a 
coarse-graining transformation W for the deformed toric 
code Hamiltonian -ff^c that explicitly preserves its local 
Z2 symmetry. 



V. COARSE-GRAINING OF A LATTICE 
MODEL WITH A LOCAL SYMMETRY 

In this section we propose a coarse-graining transfor- 
mation for lattice models with a local symm etry . For con- 
creteness, we consid er t he toric code modet^SJ deformed 
with a magnetic fielcP^as reviewed in Sect. |Hj which has 
a local Z2 symmetry, although it is possible to generalize 
the present construction to quantum double model*^ 
with a generic discrete symmetry group using the results 
of Refs. 11951 as we ll as to the more exotic context of 
string-net models^" using the results of Ref. 11961 

In order to motivate our construction, we start by dis- 
cussing two other possibilities, namely (i) the use of a 
bare coarse-graining transformation that simply ignores 
the presence of the local symmetry altogether (as a result 
of which the local symmetry may not be preserved due to 
numerical errors) , and (ii) the use of the (non-local) dual- 
ity between the Zi lattice gauge theory and the quantum 
Ising model, to apply a coarse-graining transformation to 
the Ising model instead. Then we present an overview 
of the strategy considered in this work, followed by a 
more detailed explanation of the two parts into which the 
coarse-graining transformation W splits: an exact trans- 
formation ITcxact and a numerical transformation iy num . 
Finally, we also discuss how to compose several layers of 
coarse-graining and how to build a variational ansatz fol- 
low energy states of the original model. 



A. Motivation 

Let us restate our goal. We would like to build an 
isometric transformation W, 

TT f : V £ -> V c > (46) 

to coarse-grain the deformed toric code model, 

H$ G = -j e j2A s -j m J2 B p- h *J2 a h ( 47 ) 

s p j 

defined on a lattice C made of L x L sites (or 2L 2 spins) , 
as reviewed in Sect. |Hl into an effective lattice £ . For 
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the case J e >• J m , h x , the entire low energy subspace of 
-ff^c is made of states that are invariant under the local 
symmetry, Eq. [3j We next discuss two simple options. 



1. Bare coarse-graining 

A first option is to proceed with a bare coarse-graining 
scheme that ignores the presence of the local Z 2 symme- 
try or, at most, merely exploits the global Z^ symmetry 
given by 



n< 



H, 



TC 



n< 



-"TC- 



(48) 



In this case, one could consider a transformation W 
made of several types of disentanglers and isomc- 
tries along the lines of the schemes used to study 
other two-dime nsional lat tice models with entanglement 
renormalizationSSHMSEIZ] Notice that the local symme- 
try, arguably a fundamental property of the model, will 
in general not be protected against numerical errors that 
may occur during a bare coarse-graining transformation. 
In addition, several attempts in this direction (see e.g. 
Fig. [2]) indicate that, while it is possible to coarse-grain 
the system without exploiting the local symmetry, the 
computational cost of applying the bare approach to the 
deformed toric code is prohibitively large in the decon- 
fined phase of the model, h x < /i£ nt ', due to the presence 
of large amounts of entanglement, as illustrated in Fig. 

m 



2. Duality transformation 

It might therefore be more convenient to exploit the 
local Z2 symmetry. The symmetry introduces local con- 
straints on the degrees of freedom of the model, implying 
that the vector space V c of the theory is redundantly 
large. It is possible to exploit these constraints to ob- 
tain an equivalent description in a smaller vector space, 
and to then build the coarse-graining transformation di- 
rectly on the reduced vector space. As reviewed in Sect. 
|III A[ when J e 3> J m ,h x , projection onto the low en- 
ergy subspace leads to the Zi lattice gauge theory, with 
Hamiltonian 



H 



LGT 



■£*p 



■£«* 



(49) 



where J m was set to 1 without loss of generality. Recall 
that this mo del can be mapped i nto th e quantum Ising 
modeP^EI]] as reviewed in Sect. IIID[ 



H 



Ising 



= -h. 



{p,p') 



£4 

p 



(50) 



where different topological sectors of H^. c correspond to 
different boundary conditions for -ffismg- This is the re- 
duced description we were looking for. 



Therefore a second option is to hrst transform the 
low energy sector of the deformed toric code H^ c (with 
J e 3> J m ,h x ) into the Ising model -Hi s i ng ; then to ap- 
ply coarse-graining techniques to determine the ground 
state phase diagram of the Ising model H ising j and, fi- 
nally, to translate the results back to the deformed toric 
code model by undoing the duality transformation. 

This is indeed a viable option. As a matter of fact, en- 
tanglement renormalization techniques have already been 
used to coarse-grain the quan tum Isin g model and study 
its ground state phase diagrarrt- 143 * 144 [ In particular, with 
the specific layout of disentanglers and isometries pro- 
posed in Ref. 11441 it was possible to study arbitrarily 
large lattices. The presence of a smaller amount of en- 
tanglement in the ground state of -Rising, as compared 
to the ground state of -ff^c explains why the coarse- 
graining strategy has a lower computational cost for the 
Ising model than for the deformed toric code model. 

This second approach can be extended to all those 
models with a local symmetry such that their low en- 
ergy subspace is dual to a simpler model. This is the 
case, for instance, of the quantum double modef^ for 
any discrete Abelian group Z n , whose low energy sector 
(if one allows only for certain type of excitations) corre- 
sponds to the Z n lattice gauge theory, which in turn is 
dual to then-state Potts model with a transverse mag- 
netic fielcp2] It seems, however, that for non-Abelian 
groups, an analogous duality transformation produces a 
highly non-local modeP^l. 

In spite of its simplicity, this second approach also has 
several drawbacks. Let us discuss two of them. 

Vacancy. — The duality transformation produces an 
Ising model with a vacancy that breaks translation invari- 
ance. While this may not be relevant in the thermody- 
namic limit, the presence of a vacancy in a finite system 
implies that the coarse-graining transformation cannot 
be homogeneous. As a result, for instance, the cost of 
simulations with the MERA increases from 0(log L) for 
a translation invariant system to 0(L 2 ) for an inhomo- 
geneous system. 

Weak non-locality. — Another, more fundamental lim- 
itation of using the duality is that this transformation is 
not local, as reflected by the presence of two non-local, 
boundary spins in the Ising model, see Fig. [9] We refer 
to the resulting model as being weakly non-local, since its 
Hamiltonian is still local in the bulk. A particular topo- 
logical sector of H^ c can be studied by just fixing these 
boundary qubits to some product state, corresponding 
to fixing the boundary conditions of -ffi s i ng , in which case 
the model is completely local. However, in a finite system 
it might be of interest to study processes that simulta- 
neously involve several topological sectors of the model. 
These processes can still be studied with the resulting 
Ising model by allowing the boundary spins of Fig. [9] to 
be entangled with the rest of the spins. However, each 
boundary spin is coupled to O(L) neighbouring spins. 
This level of non-locality is expected to produce large 
amounts of entanglement, with a subsequent increase in 
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computational costs, possibly rendering the approach too 
expensive. 



3. Beyond a global duality transformation 

Finally, there are two fundamental reasons why in this 
work we explore an alternative coarse-graining transfor- 
mation that is not based on simply mapping the low en- 
ergy Subspace Of -ff^C m ^° Rising- 

One is that our ultimate goal is to develop a systematic 
approach that can also be applied to locally symmetric 
models and lattice gauge theories with a more general 
symmetry group, such as e.g. non-Abelian groups, and 
not just to cases where a duality transformation exists to 
a simp ler model (as is the case for finite Abelian groups 
■Zn~)- In other words, we use the deformed toric code, 
arguably the simplest possible example, for illustrative 
purposes only. The coarse-graining transformation de- 
scribed below can be suitably generalized to quantum 
double models with arbitrary (possibly non-Abelian) dis- 
crete group (equivalcntly, to lattice gauge theories with 
arbitrary discrete gauge group) by using the results of 
Ref. 11951 as well as to string-net modela^M by using the 
results of Ref. H96l 

The second reason is that we are interested in a scheme 
that can be generalized to models with a Hamiltonian 
where the local symmetry is explicitly broken (for the de- 
formed toric code, by adding e.g. an additional magnetic 
field h y V ■ cj to -ff^c) but sucn that the l° ca l symmetry 
is still recovered at low energies. In this broader context, 
the map to the Ising model is strongly non-local, in that 
each (j y is transformed into a string operator, whereas 
the coarse-graining scheme presented here can be suit- 
ably extended in a way that locality is strictly preserved, 
as discussed in subsequent work. 



B. The strategy 

Next we describe the coarse-graining transformation 
proposed in this work. The original lattice C is trans- 
formed into an effective lattice £ , where each plaquette 
of £ is obtained by coarse-graining a block of four pla- 
quettes of C, see Fig. [18} Thus, if the original lattice £ 
is made of L x L sites, the effective lattice £ is made 
of L' x L' sites with L' = L/2. The coarse-graining is 
implemented by an isometric transformation W , 



c 



IT f :V/ 



(51) 



that maps the original Hamiltonian -ff^c hito an effective 
Hamiltonian H' 



h' = w f m r w. 



TC 1 



(52) 



What makes transformation W special is that the effec- 
tive H' exactly retains the local Zi symmetry of H^ c . 
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FIG. 18. Upper part: The coarse-graining transformation W 
transforms the lattice C into the effective lattice £ . Lattice 
jC has spins on its edges. Lattice £ also has spins on its 
edges but, in addition, has an effective free spins sitting in 
the interior of each plaquette. Lower part: Transformation 
W = VKexact VKnum breaks into an exact transformation Wexact, 
which produces an intermediate lattice C with constrained 
and free spins, and a numerical transformation if„ um , which 
coarse-grains the free spins of lattice £ 



This is accomplished by decomposing W into two parts, 



w = w. 



w n 



exact v * nunii 



(53) 



as illustrated in Figs. [T8| and [T9] The first part, Wexact, 
consists of a fixed sequence of CNOT gates (see Sect. 
Ill C ) and single-spin projections. The second part is an 
isometric transformation ITnum, to be determined numer- 
ically. 

In order to explain the role of Wexact and W, lum , it 
is convenient to distinguish between two types of spins, 
depending on how the local Z 2 symmetry acts on them. 
We say that a spin is constrained by the local symmetry, 
or just 'constrained', if it is included in the support of at 
least one star operator A s . In other words constrained 
spins are transformed non-trivially by the local symme- 
try. We say that a spin is '/ree' (that is, not constrained 
by the local symmetry) if no star operator A s acts on it. 

Notice that all the spins in C are initially constrained, 
since each of them belongs to the support of two star 
operators A s . In contrast, the coarse-grained lattice CJ 
will turn out to contain the above two types of degrees 
of freedom: constrained spins, sitting on the edges of £ , 
and free spins, located inside the plaquettes of £ , see 
Figs. [18] and 111 
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Finally, Wcxact only depends on a few structural as- 
pects of the spin model, such as the fact that £ is a 
square lattice and how the local Z 2 symmetry acts on its 
spins. In particular, Wexact is independent of the value 
of the magnetic field h x in H^ c . 

We note that Wexact is inspired in a similar transfor- 
mation proposed in Ref. 11951 for the (underformed) toric 
code, h x — 0, and that was used to show that the model 
is a fixed point of the RG flow. Here we will use W e xact as 
part of a coarse-graining strategy valid for an arbitrary 
value of h T . 



2. Numerical, approximate coarse- graining 



FIG. 19. The proposed coarse-graining transformation W 
breaks into two pieces Wexact and Wnum- W transforms the 
(constrained) spins of C into the two types of spins present 
in the effective lattice C! , namely constrained spins, repre- 
sented by triple lines, and free spins, represented by single 
lines. Wexact acts on constrained spins, whereas W nU m coarse- 
grains the free spins and acts on constrained spins diagonally 
on the a x basis. 



1. Analytic, exact coarse-graining 

Wexact is applied only on constrained spins and its role 
is to transform some of these spins into free spins, while 
leaving other spins still constrained. It maps the orig inal 



lattice C into an intermediate lattice £, Figs. 18 and 
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(Wexact 

where the vector space 
spins and free spins, 



(54) 



1 x contains both constrained 



rconst 
I 



rfrcc 
t ' 



(55) 



The details of how W OX act manages to free some of the 
spins will be explained in Sect. |VC| For now, we simply 
recall that the original lattice C contains 2L 2 constrained 
spins subject to L 2 star operators A s , Eq. [3] (of which 
1? — 1 are linearly independent, Eq. p3j). That is, there 
are only half as many star operators A s as constrained 
spins. It is therefore plausible that, by properly reor- 
ganizing the vector space Yc , transformation W C xact can 
turn some of the constrained spins into free spins. W C xact 
is also in charge of eliminating some of the spins of lat- 
tice £, namely those that after the sequence of CNOT 
gates are forced by the local symmetry to be in a fixed, 
unentangled state |+). 

An important feature of the transformation W C xact is 
that it is exact: not only it is specified analytically (as 
opposed to numerically), but it can also be exactly re- 
versed, implying that no approximation errors are intro- 
duced while squeezing the low energy, locally symmetric 
sector of -fffc m ^° ^ ne smaller vector space 
t. 



£ of lattice 



In contrast, transformation W num is mostly concerned 
with free spins. Its goal is to coarse-grain these free spins 
into (effective) free spins, while acting only in a restricted 
way (to be explained below) on the constrained spins. It 
transforms the intermediate lattice C into the effective 
lattice £', see Figs. [l8]and[l9j 



(W num )t:V £ -^V £ /, 



(56) 



where the vector space Vc also contains both con- 
strained spins and free spins, 



Y C ' = 



(57) 



W num is determined numerically, through some opti- 
mization procedure, and it is approximate, in that the 
compression of the low energy sector of H^, c into C may 
include errors, whose size and implications have to be 
monitored. 

Finally, W num depends on the specific details of Hamil- 
tonian H^ c , namely on h x . 



3. Key properties 

Before moving to a more detailed description of trans- 
formations Wexact and W num in Sects. |VC| and |VF| 
we re-emphasize the two properties that the composite 
coarse-graining map W is designed to fulfill: 

Preservation of locality. — W transforms local opera- 
tors in C into local operators in £' . In particular, we will 
see that the effective Hamiltonian H' remains local. 

Preservation of the local Z 2 symmetry. — The effective 
Hamiltonian H' retains the local Z 2 symmetry of the 
original Hamiltonian H^ c . The local symmetry is pre- 
served exactly (in spite of the fact that W contains an 
approximate, numerical part W num ) and exploited to ob- 
tain a significant reduction in computational costs (when 
compared to a bare coarse-graining strategy). 

It is the combination of these two properties, namely 
simultaneous preservation of locality and of the local Z 2 
symmetry, that distinguishes the present approach from 
the two other options mentioned earlier in Sect. V A 
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FIG. 20. Exact coarse-graining transformation W^xact- All 
spins of the original lattice £ are equivalent. However, as a 
guide to the eye, larger circles are used to denote those spins 
that will remain constrained spins in £. Left: a region of 
the original lattice £ is mapped by H4xact into a region of 
the intermediate lattice £. Blocks of four plaquettes of £ 
are mapped into single square plaquettes of £, which contain 
three free spins and a vacancy in their interior. Right: de- 
tailed sequence of the CNOT gates, denoted by red arrows, 
applied to each block of four plaquettes in order to modify 
the support of the star operators A s . For clarity, dashed ar- 
rows also show two CNOTs corresponding to neighbouring 
unit cells. At the end of the sequence, there are three types 
of spins: free spins (green), unentangled spins (in state |+)) 
and constrained spins. Unentangled spins are removed by 
single-spin projections, and are therefore not present in £. 



C. Exact transformation We Xa ct 

Transformation Wexact consists of the sequence of 
CNOT gates and single-spin projections specified in Fig. 
[20| which maps blocks of four plaquettes of C into single 
plaquettes of the intermediate lattice L. 



1. Free spins, unentangled spins and constrained spins 

The CNOT gates are applied according to a spatially 
periodic pattern that has a block of four plaquettes of 
£ as a unit cell, see right side of Fig. [20] Their goal 
is to modify the support of star operators A s . Some of 
these star operators, initially acting on four spins, end 
up acting as the operator a x on a single spin, while some 
other star operators end up acting on four spins of £. 
This produces three types of spins: 



(i) Constrained spins that belong to the support of a 
star operator A s i, where each star operator A s > acts on 
the four spins surrounding site s 1 £ C. Constrained spins 
sit at the boundaries of the square plaquettes of £,. 

(ii) Free spins, on which no star operator acts. They 
are depicted as green filled circles in Fig 



20 



Free spins 

sit in the interior of a square plaquette of C. Each free 
spin also belongs to a single-spin plaquette that emerges 
from one of the corners of a larger, square plaquette. 

(iii) Unentangled spins, constrained to be an eigenstate 
of a x . They are produced by the last three CNOT gates 
on the right side of Fig. [20]and can be identified as those 
spins sitting on an open edge (each unentangled spin is 
depicted next to a ket '|+) ). Notice that any such spin 
must be in a product state. Indeed, if the collective state 
|£) of many spins is an eigenstate of ctJ (with positive 
eigenvalue +1) for some specific spin j, then that spin j 
cannot be entangled with the rest of the spins, but must 
instead be in the state |+) (see Eq. 29), 



<I0 - 10 



10 



10 



(58) 



Notice that this property is common to all the ground 
states of the Hamiltonians H^ c independently of the 
strength of the magnetic field h x since the magnetic 
field commutes with the star operators. Transformation 
W^exact also includes projecting out any such spin, an op- 
eration that exactly preserve the state |£ re st) of the rest 
of the spins, 



10 



projection 



(+j\Z) = lOest). 



(59) 



Therefore these unentangled spins do not appear in the 
intermediate lattice £. 

In summary, constrained spins of the initial £ are ei- 
ther mapped into constrained spins of C (sitting at the 
boundary of a square plaquette), or free spins of £ (sit- 
ting at the interior of square plaquette), or they are re- 
moved. This occurs in the following proportions: 



constrained 
spins <G £ 



2 spins € £ const 

3 Spins € £free 

3 spins removed 



(60) 



Lattice L can be thought of as being made of two square 
sublattices £ = £ C onst U £f rC c- The first contains con- 
strained spins in its edges, and the second contains free 
spins (or vacancies) on its sites, see Fig. 21 



2. Transformed Hamiltonian 

Under W ex&ct , the initial Hamiltonian H^ c is trans- 
formed into a new Hamiltonian 



H = (W eX act) -fffcWcxact 



(61) 
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FIG. 21. The intermediate lattice C — £ C onst U Arcc de- 
composes as the union of a square sublattice £ CO nst with a 
constrained spin on each of its edges and a square sublattice 
£f rcc with either a free spin or a vacancies on each of its sites. 
[Shading in Cf ICC is used to indicate position relative to the 
plaquettes in £ CO nst]. 



that includes several types of terms, 



H 



s'ec 



A s 




h x K\, 



(62) 
All these terms are represented in Fig. [22] 

First, H has terms coming from the star operators A s 
of i?xc As J us t mentioned, some star operators A s are 
simply eliminated (after being transformed into a single- 
spin operator cr x ), whereas others become four-spin star 
operators A s > acting on the sites of £. Specifically, from 
every four star operators A s acting on £, only one be- 
comes a star operator A s > acting on C. 

Hamiltonian H also has terms that originate in the 
plaquette operators B p of H^ c . For every four plaquette 
operators acting on C, three are transformed into single- 
site operators a z that act on the free spins of C, whereas 
the fourth one becomes a seven-spin plaquette operator 
B n ', 



B n 



n 



(63) 



where the product includes the four spins at the bound- 
ary of plaquette p' € C and the three free spins in its 
interior. 

Finally, H also contains terms that result from trans- 
forming the magnetic field — h x ^2-a x in H^ c . Under 
Wexactj the magnetic field becomes an interaction be- 
tween free spins. There are two types of interactions: 
interactions within a square plaquette of £, and inter- 
actions across the boundary between two plaquettes of 




# ' 



FIG. 22. Depiction of the terms included in the trans- 
formed Hamiltonian H. From the structure of lattice C, 
one can readily read the new star and plaquette operators 
of H (as explained in Fig. pfh)). In addition, the magnetic 
field —h x Y^, &j in H^ c gives rise to other interacting terms 
— h x K- lrl t, which are depicted as yellow shades, (i) Two-spin 
interaction —h x a x ® a x between nearest neighbour free spins 
contained in the same plaquette. Notice that in the pres- 
ence of a vacancy, some of these terms become a magnetic 
field —h x a x acting on a free spin, (ii) Three-spin interac- 
tion —h x a x ®a x ®a x between two free spins on neighbouring 
plaquettes and a constrained spin sitting at the boundary be- 
tween those plaquettes. Again, the presence of a vacancy re- 
duces the support of some of the interaction operators, which 
end up coupling one free spin and one constrained spin ac- 
cording to ~h x o x ® a x . 



£, represented by yellow bulbs in Figs. 22 'i) and 22 ii) 
respectively. Interactions within the plaquette are of 
the form —h x a x <E) <J X and couple pairs of first neigh- 
bour free spins. [Where the first neighbour of a free spin 
is a vacancy, the interaction is reduced to a magnetic 
field —h x a x acting on the free spin.] Interactions be- 
tween free spins that belong to neighbouring plaquettes 
of C are 'mediated' by the constrained spin sitting at the 
boundary between the plaquettes, and it is of the form 
—h x a x (Si cr x <8> o x ■ [Again, where the interaction would 
involve a vacancy, it is reduced to the form —h x a x (g> a x 
acting on a free spin and a constrained spin.] 



The pattern of the interaction is suggestive. The free 
spins indeed are subject to interactions characteristic of a 
magnetically charged matter field. In this case the matter 
field transforms trivially under the action of the gauge 
group, since it belongs to Vlgt on which by definition the 
constraints act as the identity operator (see Sect 



in i 



It would then be tempting to interpret the free spins 
as an emergent matter field minimally coupled to the 
Hodge dual of the constrained spins, that play the role 
of the magnetic photons. In this way the mechanism of 
confinement as dual superconductivity of the vacuum for 
the Zi lattice gauge theory becomes manifest! 203 * 204 !. 
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D. Local implementation of the duality 
transformation 



So far we have introduced transformation W cxac t as 
a sequence of CNOT gates and single-spin projections, 
and have described the Hamiltonian H resulting from 
transforming the deformed toric code Hamiltonian H^ c . 
Let us now take a minute to analyse these results. 

Wexact has an interesting interpretation in relation to 
the duality transformation between the Z-i lattice gauge 
theory and the quantum Ising model reviewed in Sect. 
Indeed, looking back at Fig. 22 one can see that 



E3 

the three free spins inside a square plaquette of lattice C 
form a small, 2x2 quantum Ising model on their own 
(with a vacancy), in that they interact according to the 
Hamiltonian 



-h, 



E"5 



E*J 



(64) 



[where the interactions involving a free spin and the va- 
cancy are reduced to a term —h x a x acting on the spin.] 
Therefore, W exact can be thought of as implementing a 
duality transformation inside each block of four plaque- 
ttes of C. This visualization is useful. It tells us that the 
mechanism used by W eX act m order to remove degrees 
of freedom that are determined by the local symmetry 
(namely to transform some star operators A s into opera- 
tors u x acting on a single spin, which must consequently 
be in state |+), and to project out such spins) is the same 
that is used to transform the Z-2. lattice gauge theory into 
the Ising model. However, while the duality transforma- 
tion is usually applied once globally on the whole system, 
W^exact applies an independent duality transformation on 
each block of four plaquettes of C. This makes W C xact 
a fully local transformation, namely one that maps any 
local operator on C into a local operator in C 

The resulting tiny (3-spin) quantum Ising models in- 
side the square plaquettes of £, represented in Fig. 



22 



(i). 



are not independent, but rather patched together by fur- 
ther interactions represented in Fig. 
actions, 



22 li). These inter- 



(65) 



can be interpreted as Ising interactions between two free 
spins (one on each plaquette) that are mediated by the 
constrained spin sitting at the boundary between the two 
plaquettes. Specifically, depending on whether the con- 
strained spin is in state |+) or |— ), the two free spins 
interact according to —h x a x <g> a x or h x a x (g) a x . Again, 
this is reminiscent of the way the two non-local, boundary 
spin control the boundary conditions of the Ising model 
resulting from applying a global duality transformation, 
see Fig. [9j 



E. Disentangling power of Wexact 

Another relevant feature of W exac t is that it removes 
short-range entanglement present in the low energy states 
of -ffifc f° r Je ^ Jm,h x , and which has its origin in 
the presence of the local symmetry. To illustrate this 
point, we analyse the ground state of the transformed 
Hamiltonian H in the limits of small and large magnetic 
field h x , where an analytic treatment is possible. Having 
established earlier that W cxact implements some sort of 
duality transformation to the quantum Ising model at a 
local level, it is easy to anticipate that, in these two limits, 
the free spins of C will be in an unentangled state, since 
this is the case in the analogous limits of the quantum 
Ising model. 



Small magnetic field 



For h x = 0, the transformed Hamiltonian H of Eq. 
becomes 
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E^ 

s'£C 



Jnrf 




(66) 



Since in this case H is obtained by transforming the un- 
deformed toric code Hamiltonian i?xc of Eq. [2] where all 
terms commute with each other, all terms in the above 
expression also commute with each other. Therefore for 
h x = the (four linearly independent) ground states of 
H are eigenstates of each individual Hamiltonian term. 
In particular, a term —J m a z acts on each free spin in £, 
which therefore must be in state |0). Thus, for h x = 
any ground state of H factorizes as 



|^ onst ) <g> |£| ee ) 



(67) 



where |^ onst ) 

strained spins whereas \£, f ~ cc ) € 
of the free spins, 



const j g an entangled state of the con- 
"'"'"' " ^% ee is a product state 



\^°) = |0)®|0)<g>---®|0). 



(68) 



This factorization is not surprising -after all, the limit 
of small magnetic field h x corresponds to the spin polar- 
ized state of the quantum Ising model in the limit of a 
large transverse magnetic field (in the z direction). And 
yet, it is a result that has important implications for the 
purposes of this paper. It means that Woxact has trans- 
formed a robustly entangled state (the ground state of 
the underformed toric code) into a state where a fraction 
of the spins are not entangled at all. This is the basis for 
a significant reduction in computational costs. 

Following Ref. 11951 we further notice that the state 
|£const) for the constrained spins is again the ground state 
of the undeformed toric code. Indeed, after projecting 
out the free spins in state |£f rC c), the Hamiltonian in Eq. 
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66 becomes the toric code Hamiltonian of Eq. [2] defined 
now over the constrained spins of lattice £. In other 
words, the models before and after the coarse-graining 
transformation are locally identical, and the undeformed 
toric code is seen to be a fixed-point of the RG flow^l. 
For a finite h x ^ the free spins will no longer be in a 
product state, but rather entangled between themselves 
and the constrained spins. However, for small values of 
h x , the free spins are expected to still be only weakly 
entangled, an expectation confirmed by numerical simu- 
lations, see Fig. [T] Thus, one of the merits of transforma- 
tion Woxact is that it reduces ground state entanglement 
throughout the deconfined phase of the deformed toric 
mode model, see Sect. In) 



2. Large magnetic field 

We have just seen that transformation W oxact is capa- 
ble of removing entanglement from the ground state of 
H^c m the deconfined phase, a property that the present 
coarse-graining scheme will exploit to obtain a simpli- 
fied description. It would be unfortunate, however, if 
Wexact would introduce additional entanglement in the 
other phase of H^ c , namely the spin polarized phase, 
see Sect. |TT] This does not seem to be the case, see Fig. 
[TJ In the limit of large magnetic field, this is easy to see. 

Indeed, for very large h x , and using that J e ^$> J rn , h x , 
the deformed toric code Hamiltonian H^ c can be written 
as 



J, 



E 



A s - k x > f> : 



: £< 



(69) 



and its (unique) ground state is the spin polarized state 
|+) ® |+) ® • • • ® |+). It is clear from Eq. [3l] that 
the sequence of CNOT gates leave this state invariant. 
Therefore, for large h x , the ground state of H is also an 
unentangled, spin polarized state. 



F. Numerical transformation Wnum 

For a finite value of h x , the ground state of H will 
be such that the free spins are entangled with the con- 
strained spins of £. The numerical transformation W num 
coarse-grains the three free spins inside each square pla- 
quette of the intermediate lattice £ into an effective free 



spin of lattice £', see Figs. 18 and 19 

Transformation W num is required to fulfill two condi- 
tions. On the one hand it must map local operators in £ 
into local operators in £' . This can be accomplished, for 
instance, by building W, lum as a product of disentanglers 
and isometries, as explained in Sect. IV On the other 
hand, Wnum must commute with the local Z2 symmetry. 
This is achieved by requiring that, at most, it acts on the 
constrained spins by means of operators that are diagonal 
in the a x basis. For instance, a unitary transformation 

exp(-iipa x ®<j x ® a y ), ip G [0, 2tt) (70) 




FIG. 23. This simplified choice of the numerical transforma- 
tion Wnum consist of a tensor product of isometries w that 
coarse-grain the three free spins inside a squa re plaquette of 
£ into an effective free spin of £' , see Eq. 



where the first and third Pauli matrices act on free spins 
and the middle one on a constrained spin, could be used 
to remove part of the entanglement introduced by the 
coupling of Eq. |65j Since this unitary transformation 
acts on the constrained spin diagonally on the a x basis, 
it commutes with any star operator A s i acting on that 
spin, and therefore does not alter the action of the local 
Zi symmetry on £. 

In this work we use, for illustrative purposes, a simpli- 
fied transformation Wnum that is made only of isometries, 
see Fig. [23j That is, we do not incorporate disentanglers 
in the numerical part of the coarse-graining transforma- 
tion. This particular choice is clearly a limitation: lo- 
cal entanglement involving free spins is not removed by 
W num and, as a result, it will accumulate over successive 
applications of the coarse-graining transformation W, in 
a way similar to wha t is observed usin g a tree tensor 
network (TTN) ansatsPHl. i n Sect. |vj]we will see that 
this choice of W n u m limits the system sizes that can be 
addressed. 



On the other hand, restricting our attention to a nu- 
merical transformation W nU m that is made only of isome- 
tries also has several advantages from a pedagogical point 
of view. Its simplicity allows us to explicitly keep track 
of how the different terms in the initial Hamiltonian 
H^c are transformed. In addition, the absence of nu- 
merical disentanglers in W n um (whose effect has already 
been studied in other two-dimensional systems ^ 43 * 144 * 147 ]) 
makes the role of W oxac t in the whole coarse-graining, and 
in particular its disentangling power, more transparent. 
We refer to Ref. 12051 for an example of a more complex 

•'mim- 
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which again contains a number of terms, see Fig. 25 
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FIG. 24. The effective lattice £ — -Cconst U£f roe decomposes 
as the union of a square sublattice C' coust with a constrained 
spin on each of its edges and a square sublattice £f roo with 
an effective free spin on each of its sites. [Shading in £f ree is 
used to indicate position relative to the plaquettes in £const] • 



1. Effective free spins 

Let us then see how our simplified choice of W num , 
Fig. [23] transforms the system. Each isometry w in 
Wn Um simply maps the three free spins in the interior of 
a plaquette of C into an effective free spin, 



,t 



03 



'X'> 



(71) 



where C2 is the vector space of a free spin and C x > is the 
^'-dimensional space of an effective free spin. 

The dimension x' controls the degree of approximation 
introduced during coarse-graining. Ideally, it should be 
chosen large enough so that C x ' can accommodate the 
support of the reduced density matrix of the three free 
spins (see e.g. Rcf. 190). Let x' p be the dimension of this 
support. The computational cost of the approach grows 
as a power of x' ■ If; as a means to reduce the computa- 
tional cost, x' is chosen to be smaller than x' p , then the 
coarse-graining transformation becomes approximate. 

For h x — and h x = 00 the free spins are in a product 
state and x' can be chosen to be just 1, and the coarse- 
graining transformation becomes trivial. For other values 
of h x , the free spins are entangled and a larger value of 
x' must be used. 



2. Effective Hamiltonian 

Under our simplified choice of W num , Hamiltonian H 
is mapped into the effective Hamiltonian 



H' = (W num yHW num = VftHfoW, 



H' 
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(73) 



In this expression operators E z for an effective free spin 
are obtained by coarse-graining the sum of the three er 2 's 
acting on three free spins, 



w t (of + CT2 + °3:) w - 



(74) 



The operator T, x for an effective free spin is obtained by 
coarse graining the sum of the interactions inside a block, 



X^ X „i,f I ~ X ~ X I ~ X ~r X I ~r X I ~r X \ ^,, 

Zj = W ] yCF^ (T 2 + "l (Tg + CT 2 + ^3 ) w - 



(75) 



In both equations [74] and [75] the numbering of the o x 
and <7 2 operators refers to the position of the free spin as 



reported in Fig. 21 The five-spin plaquette operator B p i 
is now the product of four ct z 's corresponding to the con- 
strained spins in plaquette p' € C and of E z correspond- 
ing to the effective free spin inside this plaquette. The 
term —h x K' int collects interactions between effective free 
spins located in neighbouring plaquettes, which are con- 
trolled by a constrained spin. What is significant about 
H' is that all operators are local (their support spans a 
small number of neighbouring spins). We will see that 
this property is robust under successive coarse-graining 
transformations. 



G. Exact preservation of the local symmetry 

Let us discuss how the coarse-graining transformation 
W preserves the local Z2 symmetry. 

The local Zi symmetry acts on the original lattice C 
by means of the star operators A s , see Sect. II C As we 



(72) 



see Sect. 

have seen above, W maps one fourth of these operators 
into star operators A s i acting on the effective lattice £', 
whereas the remaining three fourth of star operators A s 
have been mapped into operators o x that act on spins 
that have been excluded from C . The roles played by 
Inexact an d W^num m this process are markedly different. 
Wexact is in charge of transforming the star operators and 
removing three fourth of them from the effective model. 
This is done analytically, so that no numerical errors are 
introduced, and with a transformation that is indepen- 
dent of the magnetic field h x , so that the way in which 
the local Z2 symmetry acts on £' is always the same. 
In contrast, Wn Um coarse-grains the system numerically 
while avoiding to modify the way in which the local Z 2 
symmetry acts. This is achieved by considering a trans- 
formation that commutes with all the star operators A s i , 
and it implies that the local Z 2 symmetry will not be af- 
fected by the possible (numerical and truncation) errors 
that Wn Um may introduce. 
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FIG. 25. 

H\ Eqs. 



Terms tha t appear in the effective Hamiltonian 
(i) Star operator A s i involving four 



72 



and 



73 



constrained spins. (IT) Plaquette operator B p r involving four 
constrained spins and an effective free spin, (iii) Operators 
£ z , Eq. [72] and E x , Eq. [75] acting on an effective free spin, 
(iv) Interactions between effective free spins on neighbouring 
plaquettes of £' , controlled by the a constrained spin. 



H. Further coarse-graining 
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FIG. 26. The coarse-graining transformation W' transforms 
a block of four plaquettes of the effective lattice £' into a 
single plaquette of a new effective lattice £" . Notice that 
lattice £" has the same local composition of constrained and 
free spins as £' . Transformation W' again breaks into an 
exact transformation W^ct , which produces an intermediate 
lattice £', and a numerical transformation Wnum- 



We can now apply a second coarse-graining transfor- 
mation 



(W'y : V/y -> V/y. 



(76) 



which again breaks into an exact part W^ act and a nu- 
merical part W^ um , see Figs. 26 and 27 



exact num' 



(77) 



Transformation W' is very similar to transformation 
W, but it acts on a lattice made of both constrained 
and free spins, whereas W acts on a lattice made only 
of constrained spins. The exact transformation W^ xact 
acts on the constrained spins of £ , and it does so by 
applying exactly the same sequence of CNOT gates and 
single-spin projections as Wexact while simply ignoring 
the free spins, see Fig. [28] to produce both constrained 
spins and fresh free spins. Then W^ um coarse-grains the 
free spins of £' and the freshly produced free spins while 
acting on constrained spins diagonally in the a x basis. 

The simplified numerical transformation W^ um consid- 
ered in this work is made of isomctries w' , 



K) T :(C x <r®(C 2 ) 
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(78) 



that map the free spins in the interior of a plaquette of 
£' (namely four effective free spins and three fresh free 
spins) into an effective free spin with vector dimension 



C x " , where again x" is the refinement parameter of the 
transformation. 

W' maps local operators into local operators in very 
similar way as W does. In particular, the effective Hamil- 
tonian H" , 



H" = (W'YH'W, 



(79) 



can be seen to contain terms analogous to those of H' . 
[To illustrate this point, Fig. 29 shows the Hamiltonian 
terms that are obtained in transforming —h x K[ nt in H' 
according to Wg Xact , and which will produce — h x K" nt .} 
With these observations, we can conclude that W also 
preserves the locality of operators and the local Z 2 sym- 
metry, and that both properties will also be preserved 
under analogous subsequent coarse-graining transforma- 
tions. 

Notice in Fig. 



26 that lattices £' and £" display 



an identical local pattern of constrained and free spins. 
Indeed, they both can be decomposed as the union of 
a sublattice for constrained spins and a square sublat- 
tice for free spins, see Fig. [24j Therefore all following 
coarse graining transformations W" ,W" , • • • have the 
same structure of W. 

This allow us to consider a larger sequence of coarse- 
graining transformations {W, W', W", ■ ■ ■ }. Each coarse- 
graining transformation halves the linear size of the lat- 
tice, reducing the number of constrained spins to one 
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FIG. 27. The coarse-graining transformation W' is com- 
posed of an exact transformation W CX act that acts only on 
constrained spins and a numerical transformation W nU m that 
coarse-graines the fresh free spins obtained from W GX act to- 
gether with the free spins in £' to produce the free spins of 
£" while acting on the constrained spins diagonally in the a x 
basis. 
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FIG. 28. First round of the CNOT gates that constitute 
Wcxact- Notice that the free spins are simply ignored by the 
CNOT gates. [In the first step, we have repositioned the free 
spins for notational convenience] . The rest of the sequence of 
CNOT gates proceeds as in Fig. [20] 



fourth and, starting with W, reducing the number of 
free spins also by one forth. Suppose that the initial 
lattice has L x L sites (or 21? constrained spins) with 
L = 2 K . After (log 2 L) — 1 = K — 1 layers of coarse- 
graining, producing a sequence {£, £'£", ■ ■ ■ £ top }, the 
original lattice has been reduced to a small lattice £ top , 
referred to as top lattice, with 2x2 sites or, equivalently, 
eight constrained spins and four free spins, see Fig. [30] 



I. Top coarse-graining transformation VF top and 
topological degrees of freedom 

The top lattice £ top is mapped into a set of just three 
spins, collectively denoted C* , by a coarse-graining trans- 



formation W top , see Figs. 30 and 31 

(W/ to P)t • 



' r to P 



(80) 



that decomposes as usual into exact and numerical parts, 
wtop _ Tj/ to P wtop 

r r r r exact r r num ' 

Transformation W top needs to be described explicitly, 
as it differs from previous transformations in a few as- 
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FIG. 29. Examples of interactions contained in the 
term (Wcxact) 1 (-hxK[ at ) Wcxact of the Hamiltonian 
( W^ct) 1 ^' Wcxact ■ (i) The interaction between two ef- 
fective free spins is mediated by two fresh spins, on which 
it acts as a x ® a x . (ii) Analogous interaction across the 
boundary of a square plaquette, which involves a constrained 
spin, (iii)-(iv) In the presence of a vacancy, the interaction 
simply acts on one spin less. 




FIG. 30. Graphical representation of £ op , the intermedi- 
ate lattice £ op , and the final lattice £* made of just two 
topological spins and a free spin. Notice the toric boundary 
conditions. 



pects. Let us start with W c °^ ct . The unit cell of four 
plaquettes on which the sequence of CNOT gates in Fig. 
l20lwas defined amounts now to the whole lattice £ top . As 



shown in Fig. 32 this implies that some CNOT gates are 
applied twice on the same control and target spins. Since 
a CNOT gate squares to the identity, such gates do not 
need to be applied. Another important difference is that 
in this case the sequence of CNOT gates does not pro- 
duce any constrained spins. Indeed, the two spins that 
would otherwise be constrained (denoted as 1 and 2 in 
Fig. 32 ) are now free. This is due to the fact that toric 
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FIG. 31. Transformation W e °^ ct produces two topological 
spins and three fresh new free spins. Transformation W^JJ, 
coarse-grains seven free spins into a single spin, while acting 
diagonally in a a x basis on the two topological spins. Topolog- 
ical sector (v\,Vi) can be specified by setting the topological 
spins in the state \vi) ® \v2)- Index a labels an orthonormal 
basis | a) G C x of states of the final free spin. Each choice 
of topological sector («i,V2) and of state \a) labels one of 
4x orthonormal states \^? v . V2 \) G Yc of the original lattice 
C, which can be reconstructed by undoing the sequence of 
coarse-graining transformations. 



boundary conditions imply that a star operator A s act- 
ing on these spins acts twice on each of them, and since 
the square of operator (a x ) is the identity operator, the 
symmetry acts trivially. These two spins are very spe- 
cial. One can check that the string operators X\ and X 2 
of Eq. [6] acting on C are mapped, under coarse-graining, 
into operators erf and o~ 2 acting on these two spins. In 
other words, the state of these spins determines in which 
topological sector the model is. For this reason, instead 
of calling these spins free, we refer to them as topolog- 
ical spins. Incidentally, on a state |£) invariant under 
the local Z 2 symmetry, Eq. |3j operators X\ and X 2 are 
equivalent to any operator obtained by multiplying them 
by star operators A s . For instance, A S X\ is equivalent 
toXi, 



(A s X 1 )\0 = X 1 A s \0=X 1 \0, 



(81) 



where we have used that X% commutes with A s and Eq. 
[3] Thus, by multiplication by star operators A s allows 
us to locally deform the original support of X\ and X 2 
(namely the non-contractible cuts C\ and c 2 in Fig. pi) 
without changing the action of these operators. In par- 
ticular, operator A S X\ above (and any such products) 
will still be coarse-grained by the present scheme into 
erf, and the same holds for deformations of the support 
of X 2 . On the other hand, the operators Z\ and Z 2 
of Eq. [8] acting on C are mapped, under the present 
coarse-graining, into operators <j\ and erf acting on these 
two spins. [Notice, however, that except in very special 
circumstances (e.g. when representing ground states of 
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FIG. 32. Lattices £ top and £ top , with toric boundary con- 
ditions, and sequence of CNOT gates and single-spin projec- 
tions included in W c °^ ct . For simplicity, the four free spins of 
C top have not been depicted. Notice that due to the boundary 
conditions, in occasions two CNOT gates (denoted in lighter 
red) are applied twice on the same control and target spins, 
in which case the do not need to be applied, since the square 
of a CNOT gate is the identity operator. 



i?TC f° r h-x — 0, which also fulfill Eq. PH) multiplying Z\ 
and Z 2 by plaquette operators B p wilTnot produce an 
operator that will still be mapped into af and erf. In 
other words, in our scheme, the support of X\ and X 2 
can be locally deformed without changing the properties 
of these operators, whereas this is not true of Z\ and Z 2l 
which must be supported on the specific non-contractible 
loops l\ and l 2 shown in Fig. [3] in order to be mapped to 
o~ z operators on the topological spins.] 

The role of transformation W^Sm i s to coarse-grain all 
free spins in C top into a single free spin described by a 
vector space of dimension x* 1 while acting on the two 
topological spins diagonally in the a x basis. For each of 
four possible choices of the topological charges (^1,^2), 
WjJuJjj is characterized by an isometry w Vl<V2 , 



(<P, 2 ) T :(C xt o P f 4 ®(C 2 r^C 



x*- 



(82) 



This completes our description of the coarse-graining 
transformations {W,W ,W" ,■ ■■ ,VF top }. By composi- 
tion, they reduce the original lattice £, made of 2L 2 con- 



strained spins with vector space (C2)' 



i,2L A 



to just two 



topological spins and one free spin with vector space 



V £ . £ C 2 



-X*- 



(83) 
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Variational ansatz 
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Let {|cc)}, a = 1, • • • ,X) denote an orthonormal basis 



'X* 



Then from states of C* of the form 



\vi) <g> \v 2 ) ® | a) e V/ 



(84) 



we can obtain a set of 4%* orthonormal states |$" „ ) <E 
V£ of the original lattice C by undoing the sequence of 
coarse-graining transformation, 



|*« )Ua ) = WW 



V7 toi >i)® |«a)(8>|a). (85) 



As it is costumary within MERA algorithms (see 
Sect. IV C), we use the set of isometric trans- 
formationFTV,^',^",--- ,W top } to define a vari- 
ational ansatz for a 4x-dimensional subspace V* of 
V£. The variational parameters, denoted a, are en- 
coded in the isometric tensors (namely disentanglers 
and isometries) that form the numerical transformations 
iW w' W" ■ ■ ■ W top \ 

\"nBrai ''numi "ram) I "numJ' 

In order to obtain an approximation of the ground 
state(s) and lowest energy excited state of Hamiltonian 
.ff^C with this ansatz, the variational parameters a are 
optimized so as to minimize the expectation value Eg of 
Hamiltonian iff c on the subspace V* C Vc, namely 



miu-EW, 



E s = tr (H* c P d ) , 



(86) 



where Pg is a projector on subspace V*. This optimiza- 
tion can be performed using the optimization techniques 
extensively described in Refs. 11451 and 11901 

The ansatz used in this work is a hybrid between a 
MERA and a TTN. The exact transformation W e xact 
contains disentanglers (sequence of CNOT gates) and 
isometries (single-spin projections) and therefore it con- 
forms to the definition of MERA. However, the simpli- 
fied choice of numerical transformation Ifn Um only has 
isometries, corresponding to a TTN. In the next section 
we refer to the ansatz as the hybrid ansatz, to be able to 
distinguish it from the proper MERA used in Ref. 
to recover the RG fixed points. 



VI. BENCHMARK RESULTS 

In this section we present the results of computations 
used to benchmark the performance of the present ap- 
proach. We have used the hybrid ansatz introduced in 
the previous section to obtain an approximation to the 
ground state(s) and several excited states of the Hamilto- 
nian iff c f° r the deformed toric code model, for lattices 
C of linear sizes L = 4, 6, 8 and 16. 

The hybrid ansatz offers an explicit representation of 
the wave function of the system, from which it is possi- 
ble to evaluate a number of quantities of interest. These 
include the expectation value of arbitrary local observ- 
ables, such as the energy, but also of non-local observ- 
ables, such as the disorder parameter (A3) and Wilson 



loops. In addition, it is possible to compute the overlap 
between different wave functions, leading to alternative 
tools to characterize the ground state phase diagram of 
the model. 

For the system sizes L = 4,6, 8, it appears that the op- 
timal choice of the numerical part of the coarse-graining 
is a simple TTN made of one layer with two isometries 
each mapping free spins into one free spin of dimension 
X top , followed by the top numerical isometry W^um- The 
refinement parameter x top i from now on denoted simply 
as Xi dominates the computational cost, which scales as 
0(x i ), see Ref. 190. Energy minimization proceeds until 
no change is observed in the first ten digits after ten opti- 
mization steps. As an example of the required computa- 
tional effort, obtaining the ground state and first excited 
state within the topological sector (+, +) for L = 4,6 
and 8 took, on a desktop computer with 2 processors at 
2Ghz, with 8Gb of Ram, a total of 1, 5 and 8 hours. 
Unless stated otherwise all the results we present here 
are obtained by fixing the value of x to x — 100. In or- 
der to determine this value, we have performed a scaling 
analysis and found that this value is enough to produce 
a consistent approximation in both the topological and 
the spin polarized phase on the 8x8 torus. The scaling 
analysis is reported in the appendix [B] 

Fig. [2] illustrates one advantage the hybrid ansatz has 
over a TTN (resulting from a bare coarse-graining trans- 
formation, see Sect. VA). In the example, given the 
same computational cost, the hybrid ansatz leads to 4 
more significant figures of accuracy for the ground state 
energy in the deconfined phase of iff c which is robustly 
entangled, see also Fig. [T] 



Local observables 



1. Low energy spectrum 



Fig. 



33 shows the low energy spectrum of iff c as a 

It includes the energy 



function of the magnetic field h a 
of the ground state and first excited state of each of the 
four topological sectors (tii,«2)> As expected, for small 
h x the ground state energies in different topological sec- 
tors are very similar, whereas they depart from each other 
for values of the magnetic fields larger than h x ss 0.3, in 
which case the ground state in the sector (+, +) becomes 
the global ground state. 

The low energy spectrum of iff c can be obtained by 
minimizing the ansatz in either the whole vector space or 
within each specific topological sector {v\,V2), which is 
achieved by fixing the topological spins to state |ui)®|f2)- 
Since the second option needs to deal with less low en- 
ergy states at a time, it is generally more economical. 
In addition, in the spin polarized phase, where ground 
states in different topological sectors have markedly dif- 
ferent energies, restricting the minimization to a single 
topological sector, e.g. (+,—), is important in order to 
get an accurate representation of its ground state. In- 
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FIG. 33. Lowest eigenvalues of H^ c as a function of the mag- 
netic field h x , for a system of linear size L = 8, corresponding 
to 8 2 x 2 = 128 spins. The two lowest energies in each topo- 
logical sector (111,1)2) are plotted. The lowest energy in each 
sector is represented with a solid line while the second lowest 
with a dashed line. Different colours refer to different sectors. 
The energies in sector (-,+) are identical to the energies in 
sector (+,-)■ The inset shows the gap between the energy of 
the ground states in sector (vi,Va) and in the sector (+,+). 
These gaps are very small for values of the magnetic field h x 
smaller than some critical value h x Tlt — 0.3285^-^ (denoted 
by a vertical line) and much larger for larger magnetic field 
h x , see Fig. |34| These results were obtained with the hybrid 
ansatz with \ — 100. 



deed, for large h x the ground state |$+,_) is no longer 
a low energy state of H^ c when the whole vector space 
is taken into account (as shown next, |4>( +i+ )) has 0(L) 
lower energy) , and will not be properly captured through 
an unrestricted energy minimization. 



2. Gaps between topological sectors 




FIG. 34. Upper panel: Energy gap AE between the ground 
states of the topological sectors (+, — ) and (+, +) (sec Eq. 
87 1 as a function of the linear size L for h x = 0.2 and 0.3, 
corresponding to the deconfined phase of the deformed toric 
code model. In accordance with Eq. |12[ AE decays exponen- 
tially with L. These results were obtained with the hybrid 
ansatz with \ — 100. Lower panel: Energy gap AE between 
the ground states of the topological sectors (+, — ) and (+, +) 
(see Eq. 871 as a function of the linear size L for h x = 0.5 



and 0.55, corresponding to the spin polarized phase of the 
deformed toric code model. In accordance with Eq. |I8| AE 
grows linearly with L. These results were obtained with the 
hybrid ansatz with \ = 100. 



Fig. [34] studies the gap 

AE = E , _ - E 



+.- 



(87) 



between the ground states of the topological sectors 
(+, — ) and (+, +) as a function of the system size. We 
obtain a clear confirmation of the exponential decay an- 
ticipated for the deconfined phase, see Eq. [T2j as well as 
of the linear growth of the gap with system size in the 
spin polarized phase, see Eq. 
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3. Gap within the topological sector (+, +) 



Recall that the ground state |<J>+,+) is the absolute 
ground state of -ff^c ^ or an arbitrary value of the mag- 



netic field h x . Fig. 35 shows the energy gap AE + ._ 



within the topological sector (+, +), as a function of h x , 
and for L = 4, 6, 8. Notice that the minimum of AE +y+ 
as a function of h x closes roughly as 1/L. This minimum 



occurs for a value of h x that diminishes with L and a 
rough estimate of the critical value h x Tlt in the thermo- 
dynamic limit can be obtained by a large L extrapolation. 
However, a much more accurate estimate of h x rlt is ob- 
tained by drawing the curves of LAE+ + for increasing 
values of L in the range L = 2,4,6,8. The position of 
intersections between those curves obtained at consecu- 
tive L (i.e L = 2 with L = 4, L = 4 with L = 6 and 
L = 6 with L = 8) produces a sequence of h x nt (L) that 
following the seminal work of Nightingale in Ref. 12061 is 
expected to converge to the location of the fixed point 
as h x nt (L) = h c x rlt - A/L 4 , see Ref. [Ml By fitting 
this expression to the sequence of intersections reported 
in the lower panel of Fig. [35] we obtain an estimate of 
the critical point ^"* = 0.3267(5). This value is off the 
exact value 0.328^^ by only 0.5%. The systematic er- 
ror induced by limiting x to \ — 100 indeed pushes the 
transition towards smaller h x since it limits the amount 
of entanglement in the ground state, however its effects 
are surprisingly small. 
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FIG. 35. Upper panel: Gap AE +t+ between the ground 
state and first excited state of the topological sector (+, +), 
as a function of the magnetic field h x , and for system sizes 
L — 4, 6 and 8. Lower panel: Product of the systems size 
L times the energy gap A^^, plotted as a function of the 
magnetic field h x (compare with Fig. |35[ ). An estimat e of 
0.3267(5) for the critical magnetic field h x rit = 0.3285(1^ 
(denoted by a solid vertical line) is obtained by ex trap olating, 
using the phenomenological renormalization grourF^, the se- 
quence of intersections among curves obtained from systems 
with consecutive sizes (2 — 4, 4 — 6 and 6 — 8) identified by dot- 
ted vertical line s in the plot using the known L dependence 
for the sequence^ 221 . 
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FIG. 36. Disorder parameter (X3) as a function of the mag- 
netic field h x , for sev eral lattice sizes. The critical magnetic 
field h x Tit = 0.3285(1^ is denoted by a vertical line. The 
inset shows a log-log plot of (X3) at h x = h c x A as a function 
of L, from which we obtain an estimate of 0.96(5) for 2j3/v, 
whose current estimates are around 1.09^2. 



f3 = 0.33 ± 0.02 for the critical exponent j3 that should 
be compared with the best estimates coming from Monte 
Carlo simulations /3 — 0.32652(15)^1. It is important 
to notice that a direct measurement of this operator in 
Monte Carlo simulations is very challenging. Only re- 
cently, in the context of U(l) LGT, an algorithm has been 
proposed to directly compute these kind of operators^!!. 



2. Wilson loops 



Non-local observables 



Disorder parameter 



Fig. 36 shows the disorder parameter (X3), introduced 
in see Sect. |IIB| which measures the formation of a con- 
densate of magnetic vortices (or magnetic monopoles), 
as a function of the magnetic field h x . Results for dif- 
ferent system sizes L = 2, 4, 6 and 8 are presented. (X3) 
can clearly be used to distinguish the deconfined phase, 
where it vanishes, from the spin polarized phase. Notice 
that the drop of (X3) to zero for diminishing magnetic 
field h x near the critical point becomes sharper with in- 
creasing system size L. For L = 8, we obtain an estimate 
h c x rit = 0.327 ± 0.05 for the critical magnetic field and 
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Fig. [37] shows the scaling of Wilson loops, see Sect, 
for the ground state of the topological sector (+, +). 
In this case, a lattice of 16 x 16 sites, or 16 2 x 2 ~ 512, 
spins was considered for a magnetic field h x — 0.2 (de- 
confined phase) and h x — 1 (spin polarized phase). For 
these values of the magnetic field the system is away from 
the critical point, as reflected in a drop in the amount of 
ground state entanglement, see Fig. [T] This allows the 
hybrid ansatz to reliably represent much larger lattices 
with just x = 100. 

For h x = 0.2 our results confirm that a Wilson loop de- 
cays exponentially fast with the size of the loop according 
to a perimeter law, as it is characteristic of the decon- 
fined phase. Instead, for h x — 1 the decay is exponential 
in the area encircled by the loop, as it is characteristic of 
the spin polarized phase. 
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FIG. 37. Upper panel: Wilson loop (Z[t\) as a function of 
the perimeter p(l) of loop I, on a 16 x 16 lattice with magnetic 
field h x = 0.1, corresponding to the deconfined phase. Notice 
the exponential decay of {Z[l]) as a functio n o f the perimeter 
p(l) of loop I, or perimeter law, see Eq. 17 Lower panel: 
Wilson loop {Z[l]} as a function of the area a(l) of loop I, on 
a 16 x 16 lattice with magnetic field h x — 1, corresponding 
to the spin polarized phase. Notice the exponential decay of 
{Z[l\) as a function of the area a(l) of loop I, or area law, see 
Eq. [201 
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their boundaries. Indeed, the ground state wave func- 
tion is somewhat similar within a phase and experiences 
a radical change when the system undergoes a phase tran- 
sition, with these two facts being captured by the ground 
state fidelity. Interestingly the characterization of phase 
boundaries using the fidelity does not require the knowl- 
edge of the order parame ters. He re we will focus on the 
intensive fidelity f(h x , h' x ^pHHUl ; defined through 



log f(h x> ti x ) 



1 
L 2 



\o S (\($(h x )\$(K))\) 



The present approach yields an explicit representation 
of the ground state wave function, from which comput- 
ing overlaps is straightforward (this is true of the hybrid 
ansatz used in this work but not in the m ore general 
case where Wn Um also contains disentanglers ?05 ). Fig. 
38 shows the logarithm of the intensive fidelity 



f o (h x ) = f(h x ,0) 



(89) 



between the ground states of -ff^c m * ne topological sec- 
tor (+, +) for h x — and another value h x > 0; as well 
as the logarithm of the intensive fidelity 



/oo(M = f(h x , oo) 



(90) 



between the ground states of H^ c in the topological sec- 
tor (+, +) for h x = oo and another value h x > 0. Fidelity 
fo(hx) shows a markedly different behaviour for small 
magnetic fields roughly smaller than h x rlt (corresponding 
to the deconfined phase, to which |$(0)) also belongs to) 
and for large magnetic fields, where the fidelity vanishes. 
An analogous behaviour is also observed for foo(h x ). 
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FIG. 38. Logarithm of the intensive fidelities fo(h x ) and 
faotha,) as a function of the magnetic field h x , for different 
values L = 4, 6 and 8 of th e system size. The critical magnetic 
field hf H = 0.3285(1^ is denoted by a vertical line. 



2. Topological fidelities 

For h x = 0, the ground states in two topological sec- 
tors only differ in the expectation value of operators X x 
and X2, see Sect. |II A| One can map e.g. the ground 
state |$+,+) of topological sector (+,+) into the ground 
state |$+,-) of topological sector (+, — ) by just applying 
operator Z 2 of Eq. [8j namely |$+,_) = Z 2 \$ + , + ). There- 
fore (4> + . + |Z 2 | < I > + ._) = 1 and the topological fidelity /+.'+ 
between these two sectors, defined as 



log(/J 



1 

I? 



log «<!>+,+ |Z 2 1 $+,_)), (91) 



Wave function fidelities 



Ground state fidelities 



The overlap or fidelity ($(/i x )|$(/i^)) between the 
ground state of H £ c for two values h x , h' x of the mag- 
netic field can be used as an alternative way to detect 
the presence of phases in the model and the location of 



is maximal, f+'Z = 1, for h x = 0. More generally, the 
topological fidelity between sectors (fi,^) and (v[,v 2 ), 
defined as 



1 



iog(/£;£) = -logi^^Kz.r {Z2) W2 1$^), 



where W; is if V; = v 1 , and 1 otherwise, also fulfils 



(92) 



f v t ,v 2 

JVl,V2 



(h x = 0). 



(93) 
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FIG. 39. Logarithm of the intensive topological fidelity fl'I 
as a function of the magnetic field h x . In the deconfined phase, 
/+'+ remains large (its logarithm remains close to zero), indi- 
cating that the ground states |$+,+) and |$+,_) differ mostly 
on their topological spins but are otherwise very similar. This 
trend changes for larger magnetic fields h x . As we enter the 
spin polarized phase, the ground states |$+,+) and |$+,-) 
are structurally different, and modifying the topological spins 
of e.g. |$_|_._) does no longer produce a good approxima- 
tion to |$+,+). Left: data collapse occurs in the deconfined 
phase when we plot Llog(/7T'7). The critical magnetic field 
h% rit = 0.3285(1^ is denoted by a vertical line. Right: in- 
stead, plotting log(/j?7) leads to data collapse in the spin 
polarized phase. This behaviour is expected as explained in 
the main text below equation |94| 



For h x 7^ Eq. [93] is no longer expected to hold. 
In particular, in the limit of large magnetic field h x a 
simple analytical characterization of the ground states 
|$ Wl ,i! 2 ) exists and can be used to show e.g. that 
|($+,+ \Z 2 |$+,_) | = L- 1 / 2 , so that /+'" vanishes for 
large L. More generally, 



f v 1 ,v 2 

Jvi,V2 







(h x — oo, L = oo). 



(94) 



For intermediate values of h x , the topological fidelity is 
expected to remain close to one in the deconfined phase 



and sharply decay in the spin polarized phase. Fig. 39 
for log(/+ ^7) as a function of h x shows that this is in- 
deed the case. In addition, the data for different values 
of L collapses into a single curve in the deconfined phase 
if log(/ + : +) is multiplied by L, whereas in the spin po- 
larized phase the data collapses into a single curve when 
log(/7'7) is directly plotted. This is reminiscent of the 
two possible behaviour of Wilson loops, namely perime- 
ter versus area law, see Fig. 37 After all, operators 



Z\ and Zi can be understood as Wilson loops with non- 
contractible support. 



VII. DISCUSSION 

In this paper we have proposed a coarse-graining 
transformation for lattice models with a local symmetry 



that simultaneously preserves locality and the symmetry, 
while exploiting the latter to significantly reduce com- 
putational costs. This coarse-graining transformation, 
made of an analytical part W^xact and a numerical part 
Wnumj gives rise to a variational ansatz for the ground 
state(s) and low energy states of the lattice model. Here 
we have focused on a simplified form of the ansatz, a 
hybrid between a TTN and the MERA, in which W„ um 
does not contain disentanglers (although Wexact does). 
The computational cost of the resulting approach grows 
exponentially with the linear size L of the lattice, severely 
restricting the size of systems that can be analysed, al- 
though it permits to consider systems well beyond the 
scope of exact diagonalization techniques. 

By adding disentanglers to Wnum, it is possible to re- 
move more entanglement from the ground state. The 
resulting ansatz is a proper MERA, with which much 
larger systems can be addressed. In addition, as ex- 
plained in Ref. 12051 under coarse-graining transforma- 
tions the ground state of H^ c for a small magnetic field 
h x can now be seen to flow back to the h x = fixed point 
of RG flow, which has a richer, local Z 2 x Z 2 symmetry, 
whereas the ground state for h x > /i"' 1 * flows to the spin 
polarized RG fixed point characterized by an infinite h x . 



A. Generalizations 

The coarse-graining transformation proposed in this 
paper can be generalized to more complex settings in a 
number of ways. 



1. Static electric charges 

Throughout the paper we have assumed that the state 
of the lattice model satisfies the star constraints of Eq. 
|3J corresponding to the absence of electric charges. How- 
ever, the coarse-graining scheme and related ansatz can 
be readily adapted to the case where one or several star 
constraints are violated on specific sites, provided that 
the state £ of the lattice is constrained to satisfy 



M0 = -\0 



(95) 



on those specific sites. Indeed, the only change that is 
required in the coarse-graining scheme is that when the 
star operator is eventually transformed by a sequence of 
CNOT gates into a single-spin operator — cr x , the spin on 
which it acts, forced now to be in state |— ) instead of 
|+), is projected out accordingly. 

In this way one could study the ground state of the 
system in the presence of two electric charges separated 
a distance /, and obtain the ground state energy as a 
function of I, from which the string tension can be eval- 
uated. 
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2. Magnetic field in the xz-plane 

If the toric code Hamiltonian is perturbed by adding a 
magnetic field — ^Ah x a x + h z a z ) on the xz-plane, a lat- 
tice model with local Z 2 symmetry can still be recovered 
by adding dummy spin variables, corresponding to mat- 
ter field, on the sites of lattice £, see Ref. 12171 In that 
case the local Z 2 symmetry can still be exploited, but 
in a trivial sense (the dummy spin variables are factored 
back out as part of the coarse-graining transformation) 
that appears to have no advantage to the case of an ar- 
bitrary perturbation, which is considered next. 



3. Arbitrary perturbation of the toric code 

In this work we have focused on a perturbation of the 
toric code Hamiltonian H^q in Eq. [2] that breaks its 
Z 2 x Z 2 local symmetry for h x = 0, generated by both 
star operators A s and plaquette operators B p , into a lo- 
cal Z 2 symmetry for h x ^ generated by star opera- 
tors only. A generic perturbation of Hamiltonian -ffxc 
will break the local symmetry completely. However, as 
recently proven by Bravyi, Hastings and Michalakis in 
Ref. 12181 topological phases are robust under arbitrary 
(sufficiently weak) perturbations of the Hamiltonian. 

For a perturbation that completely breaks the local 
Z 2 x Z 2 symmetry of the toric code model, say a weak 
magnetic field in a direction r = (r x , r y ,r z ) different from 
directions ±x and ± J 219 I 22[ M we can modify Wexact so 
that it no longer projects spins into |+) state. The se- 
quence of CNOT gates in Fig. [20j which produces six 
spins per plaquette on which either operator —J e a x or 
operator —J e o z is acting, is now followed by an extended 
numerical transformation W, lum that acts on all these six 
spins to produce an effective spin inside each plaquette of 
C . For small values of the perturbation, the sequence of 
CNOT gates will still map a robustly entangled ground 
state into a weakly entangled ground state, lowering sig- 
nificantly the computational cost of the approach, while 
the use of disentanglers as part of W nVLla is expected to 
again produce a flow back to the RG fixed point given by 
the undeformed toric code Hamiltonian H^ c . 



4- Quantum, double models and string-net models 

The analytic transformation W cxac t used in this work 
is, in essence, equivalent to the RG transformation pro- 
posed in Ref. 11951 to show that the toric code model is 
a fixed point of the RG flow. In that work, analogous 
analytical transformations were proposed also for gen- 
eralizations of the toric code, which corresponds to the 
quantum double of the Z 2 group, to models based on the 
quantum double of any discrete group Cfi^. In the ap- 
propriate regime, the low energy sector of the quantum 
double model with group G corresponds to a lattice gauge 
theory with the same gauge group. A similar analysis was 



subsequently carried forward in Ref. 11961 for string-net 
models^ 21 . 

The analytical transformations described in Refs. 11951 
and 11961 for the fixed-point Hamiltonians of quantum 
double and string-net models map robustly entangled 
spins locally into spins that are in a product state and 
can therefore be factored out (or projected out) from the 
coarse-grained system. They can again be used as the 
basis for a generalized analytical transformation W eX act, 
which needs to be supplemented with a numerical trans- 
formation Wnum when the fixed-point Hamiltonian is de- 
formed with an arbitrary perturbation. In this way, we 
obtain a coarse-graining scheme for deformed quantum 
double models (equivalent in the appropriate regime to 
lattice gauge theories) and string-net models, as well as 
corresponding tensor network ansatze for their ground 
state wave functions. 



5. Renormalization of PEPS 

The use of an analytical transformation W CX act as part 
of a coarse-graining procedure can be also exported to 
the domain of approaches based on PEPS (also referred 
to as TPS). This can occur in two different contexts. 

On the one hand, given a PEPS representation for a 
state |f) (e.g. the ground state of Hamiltonian H^ c ) 
with a local Z 2 symmetry, an important task is to evalu- 
ate the tensor network corresponding to the scalar prod- 
uct (f|f) (and related quantities), which plays a central 
role in the computation of expectation values of observ- 
ables and various optimization algorithms. The tensor 
network for (f |f) can be evaluated by a number of dif- 
ferent methods , such as corner transfe r matrix (CTM) 
methods^ 211222 ! MPS technique^- 1 26 ' 127 1 or th e tensor en- 
tanglement renormalization grourJ 129 | 132 J12al (TERG). All 
these methods proceed by coarse-graining the tensor net- 
work for the scalar product (f |f) (as opposed to coarse- 
graining a tensor network for the state |f)), with a cost 
that depends on the amount of entanglement in the sys- 
tem. This cost can again be significantly reduced if the 
coarse-graining incorporates analytic local pre-processing 
along the lines of the transformation W c ^ ac t discussed in 
this work. 

On the other hand, a PEPS representation for the same 
state |f) can also be coarse-grained as a wave function. 
This process, termed wavefunction renormalization in 
Ref. I223[ i s thoro ughly equivalent to enta ngle ment renor- 
malizatiori 141 ! 145 ^ and produces a MERA13H as a result. 
Therefore, the exact transformation Wexact discussed in 
this paper can also be used in the context of wavefunction 
renormalization. 



B. Conclusions 

Tensor network algorithms offer a variational approach 
to strongly interacting systems on a lattice that is free 
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of the sign problem and can therefore be applied to sys- 
tems that are beyond the reach of Monte Carlo sampling 
techniques, such as frustrated antiferromagnets and in- 
teracting fermions. In this paper we have explored the 
use of tensor network techniques in the context of lattice 
gauge theory and, more broadly, systems with topolog- 
ical order. We have explain how to incorporate a local 
symmetry into a coarse-graining scheme and the result- 
ing variational ansatz, for the simplest non-trivial case of 
the Z^ lattice gauge theory. 

We envisage that proper generalizations of the results 
presented here, possibly along the lines of those discussed 
in Sect. |VII A| will constitute the basis for future numeri- 
cal simulations of lattice gauge theories and more general 
models with topological order. 
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Appendix A: Duality transformation 

In this appe ndix w e review the well-known duality 
transformation^ 2 ! 154 ! between Z2 lattice gauge theory 
and the quantum Ising model in two spatial dimensions 
as described in Sect. |IIID| see FigjHJ The transformation 
is expressed in terms of CNOT gates in order to high- 
light its similarities and differences with the analytical 
transformation W CX act used in this work. 

As indicated in Figs. |40]|44| for the case of a lattice 
of 4 x 4 sites with periodic boundary conditions in both 
directions, this transformation can be implemented as a 
sequence of CNOT gates as indicated. Fig. |40| represents 
six stages (i)-(vi) of this transformation, with stages (i) 
and (vi) corresponding to the Z2 lattice gauge theory and 
the quantum Ising model respectively. 

Fig. 41 indicates the sequence of CNOT gates that 
transform stage (i) into stage (ii). This sequence pro- 
gresses from top to bottom and can be divided into four 
columns of CNOT gates, where gates in different columns 
commute. 

Fig. [42] indicates the sequence of CNOT gates that 
transform stage (ii) into stage (Hi). Again, the sequence 
progresses from top to bottom and can be divided into 
four columns of CNOTs gates, where gates in different 
columns commute. 



Fig. 43 explains how to transform stage (Hi) into stage 
(iv). Notice that stages (iv) and (v) arc equivalent, since 
they only differ in where the single-spin plaquettes are 
connected to the rest of the lattice - that is, they corre- 
spond exactly to the same Hamiltonian. 



FIG. 40. Duality transformation between the Z2 lattice gauge 
theory and the quantum Ising model, broken into six stages 
(i)-(vi). 
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FIG. 41. CNOT gates that transform stage (i) into stage (ii). 
The sequence proceeds from top to bottom. Only two of the 
three required steps are depicted. In a L x L lattice, L — 1 
such steps are required. 



Appendix B: Assessment of the precision of the 
numerical results 



Finally, Fig. 44 shows how to transform stage (v) into 



stage (vi.), corresponding to the Ising model. 



Section [Vl| presented a number of numerical results for 
the deformed toric code Hamiltonian -ff^c on a torus of 
linear size L = 4, 6 and 8, obtained with a hybrid tensor 
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FIG. 42. CNOT gates that transform stage (ii) into stage 
(Hi). The sequence proceeds from top to bottom. Only two 
of the three required steps are depicted. In a L x L lattice, 
L — 1 such steps are required. 



FIG. 44. CNOT gates that transform stage (v) into stage 
(vi). Three required steps are depicted. In a L x L lattice, 



L — 1 such steps are required. 
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FIG. 43. CNOT gates that transform stage (Hi) into stage 
(it)). Notice that stage (Hi) contains four columns of single- 
spin plaquettes, with three spins on each column. Instead, 
stage (iv) contains four columns of single-spin plaquettes with 
four spins on each column, except for the rightmost column, 
which only contains three spins and a vacancy. The sequence 
proceeds from left to right. Only two of the three required 
steps are depicted. In a L x L lattice, L — 1 such steps are 
required. 



network with a fixed value of the refinement parameter 
X = 100. Those results successfully reproduced both 
quantitatively and qualitatively the main properties of 
the ground state of the system in the deconfincd and spin 
polarized phases and of the second order phase transition 
between them. In this appendix we perform a scaling 
analysis with respect to %, to show that the particular 
choice x = 100 used in Sect. |VI| did not have a significant 
effects on the numerical results. 



We start by reminding that the hybrid tensor network 
ansatz is, as far as its numerical part W, lum is concerned, 
a tree tensor network TTN, namely one used to approx- 
imate the ground state of the lattice model after being 
transformed according to W exa ct- As discussed in Ref. 
11901 a sufficiently large value of x wn l reproduce the 
ground state to arbitrary accuracy, but this value will 
have to be exponentially large in the linear size L of the 
lattice. 

Our scaling analysis consists of two parts. First we 
study, in a torus of linear size L = 8, the convergence of 
the ground state energy as a function of 1/x, see Fig. 45 
for specific values of the magnetic field h x , namely h x = 
{0.1, 0.2} (deconfined phase), h x — {0.5, 0.6} (spin polar- 
ized phase) and h = 0.33 sa h x rlt (critical point). From 
these results, lower and upper bounds to the ground state 
energy are obtained by observing that the ground state 
energy is a monotonic function of 1/x, E g(l/ x), with 
positive, monotonically increasing derivativ e 190 !. Assum- 
ing that this trend continues all the way to 1/x = 0, this 
implies that Eq(x — 200) provides us with the best upper 
bound to the ground state energy, while a lower bound 
can be found with a linear extrapolation of Eq(x = 150) 
and E (x = 200) (the two best results) to 1/x = 0. 
For all the values of the magnetic field considered these 
bounds provide narrow windows (around 0.01% of the 
actual value). 

Once we have established the level of convergence of 
the ground state energy, we look at the dependence of 
other observables with x 
order parameter (X3 



Figs. 



46 and 47 show the dis- 
and the gap A_E( + . + ) in the (+, +) 
sector, respectively. In contrast with the ground state en- 
ergy, which was monotonic in 1/x, we see that (A3) and 
A£V_|_ + ) behave more erratically as a function of l/x- I n 
this case, we can use the variations in value as a function 
of l/x as an estimate of the error introduced by finite x- 
Importantly, this error is much smaller than the varia- 



:W 



tion of (X3) and AE/ +!+ ) as a function of the magnetic 



field h x . Further description is left to the captions of the 
figures. 
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FIG. 46. Scaling of the gap AE++ for the deformed toric 
code on the 8x8 torus as a function of 1/x for different values 
of the magnetic field in the topological phase, h x — 0.1 and 
h x = 0.2 (top panel), at the critical point , h x rlt — 0.3285 
(central panel), and in the spin polarized phase, h x = 0.5 
and h x — 0.6 (lower panel). The dependence of the gap as a 
function of 1/x is not uniform as in the case of the ground 
state energy and thus we cannot extract rigid bounds on its 
value. However due to the small variations of the gap as a 
function of \ i n the range 50 < ^ < 200 (in the worst case of 
the critical point of about 1% of its value) compared with its 
large variation as a function of the magnetic field h x , we can 
neglect the dependence on \ °f the gap provided we accept 
an uncertainty of about 1% on the actual gap value. 
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FIG. 45. Ground state energy Eq for the deformed toric code 
on a 8 x 8 torus as a function of the refinement parameter 
1/X for different values of the magnetic field, i) h x =0.1 and 
h x = 0.2 in the topological phase, ii) h c x r%t = 0.3285 at the 
critical point, iii) h x — 0.5 and h x = 0.6 in the spin polarized 
phase. Results of the simulations in the range of refinement 
parameter 50 < \ < 200 are shown with blue circles. We 
estimate upper and lower bounds for the ground state energy 
using the analysis of Ref. 11901 The upper bound indicated 
by a black horizontal line is the value of Eq at x = 200 while 
the lower bound is shown by a green line. In the plots we also 
write the exact numerical values for both upper and lower 
bounds. They provide narrow windows (in general around 
0.01% or less of the exact value) in which the exact ground 
state energy is contained. 
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